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ABSTRACT 

A program for control system modeling, simulation, and 
realization useable by the novice programmer was developed 
for interactive use on a microcomputer. The program uses two 
software packages; PL/1-60 from ligital Research and LEASE II 
from Ashton Tate. 

The paper covers a complete example for modeling and 

simulation of a minimum time response ripple free controller 
2 

for a 1/S Plant. The program uses two sample rates, one for 
modeling the computer and the other for modeling the plant. 
Variable delays due to computer computation are includea and 
corrected semi-transparent to the control designer. 

This thesis covers two different fields of interest, one 
for the experienced control system designer yet a novice 
programmer, and the other for the novice control designer yet 
an experienced computer programmer. Key programming concepts 
include using a business database to automatically program a 
control problem requiring a scientific language. 
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I. INTROEUCTICN 



A. CONCEPT 

The intent of this thesis was to develop a corr.puter 
program that an experienced controls engineer could use to 
help him model, simulate, and realize control s^zster.s with 
the convenience of a aeshtop mi cro zcmputer . Tte control 
engineer will not need to be an experienced computer 
programmer. In order to accomplish the above requirements it 
is necessary that the program be easy to use, effectively 
requiring that the program write the com.pu ter code itself 
with little help from the control designer. This is commonly 
referred to as automatic programming. In order for tre 
program to be a useable tool for the controls designer the 
system must also produce a good model of the system tiat is 
being designed. Obviously all control systems cannot he 
considered in this type of project. The anproach taA'er in 
this thesis was to use control systems that have error 
si finals designed to be equal to zero or differ from zero by a 



constant. Other guidelines of this project include: the 
equations must be simple difference equations, the system 
must allow for calculation delays within the computer and 
these delays must be corrected for automatically or at least 
allow for the possibility that the designer can correct for 
the delays, and the program must be able to model the plz.nt 



between the sample periods cf the computer which basically 
means there are two sample periods, one for the oomputer and 
one for the plant. Since the plant cannot be programmed as 
analog it must be discretized. Using a faster sample rate for 
the plant allows simpler equations yet maintains a good 
representation for the plant. Additionally, the program must 
be able to model analog to digital and digital to analog 
conversions. After the model has teen designed and the 
simulation run, the results that are available from the 
program must be easy to use and easy to interpret by the 
controls engineer. In order to meet all of these 
requirements, most of the program development itself must be 
transparent to the controls engineer. The approach taken was 
for the control engineer to provide only equations. function 
names, and function relationships by providing the names 'f 
other functions used as inputs for each equation. 

B. ASSUMPTIONS AN! H2CUIREMENT3 

The assumptions and requirements are stated below, ^irst. 
as mentioned above, the program must oe relatively easy to 
use. The program, must have high precision math capability 
since it was designed to be used with tontrcl =yste"s. Speed 
of compiling and operation was not considered critical at 
this point since micropro ressors are becoming smaller and 
faster which will solve any realization problems in the rear 
future. A major consideration is that once ‘he simulations 
have been run and the results are acceptable to the designer 
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then the simulation equations should be easily converted into 
realization equations. Basically this means the simulation 
equations in the computer should be easily converted into 
hardware equations that can be used to control a real plant. 
In order to accomplish this task it was decided that the 
program must develop source code that could be compiled by a 
high level language compiler for either computer simulation 
or hardware implementation. This approach makes the prcgram 
transportable and hardware independent. 

C. IMPLEMENTATION 

This project was implemented entirely on -microprocessor 
based systems to show that it is feasible and practical to 
use these small systems for major development. The present 
system is slightly constrained due to the length of compile 
time. All programs were developed using a Z8C micronro cessc r; 
however, the code is written in high level source languages 
therefore the programs will also run on the simpler SOSO 
based microcomputers or on the more complex sixteen bit 
microprocessors such as the 80S6. 

The languages chosen were PL/1-S0 from ligital Research 
oecause of the math and hardware capabilities ff the language 
allowing for control of analog to digital and other orocesses 
inside the computer, and DBASE II from Ashton Tate Ahi:b 
allows easy and high level manipulation of databases and file 
systems. One of the major points of this thesis is to 
manipulate functions and functional relationshios as 
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databases. Once the relationships have been developed, the 
database program writes PL/1 source code from the database 
information. The PL/1 source code is then compiled by the 
PL/1 compiler and the simulation run on the resident 
hardware . 

The present programs and the test example in this report 
provide simulation only and send the results in floating 
point format to a disk file. The floating point file can then 
be converted to display the information in tabular cr 
graphical form such as shown at the back of this report. The 
information can be displayed in a variety of forms such as on 
display consoles, printers, plotters, or graphic terminals. 
The real power for design is in the use of graphical 
equipment to provide the control engineer with the entire 
picture of the system under test. The Hewlet Packard 9872B 
plotter and 2647A graphics terminal and the Intercolor 8363 
graphics terminal were used to show that the data produced 
was basically machine independent. The code can be compiled, 
run, and displayed on a variety of different machines. 



II. GENERAL EISCUSSION 



A. CONTROL ENGINEER INFORMATION 

The control system engineer does not need to know how the 
computer programs work in order to use the system. It is 
assumed that the control engineer is a novice programmer with 
some knowledge of the CP/M operating system and limited 
background in program writing? however, it is assumed he is a 
Knowledgeable control designer. This means he is fairly well 
versed in taking control systems and producing discrete 
system equations from the continuous control equations. These 
equations are designed and optimized by whatever means are 
available either by hand or on a mainframe computer. 
(Microcomputers will scon be able to io major system 
development equations.) Once the equations are developed the 
designer enters them into a microcomputer and tests the 
simulations. If the simulations are accentable then it should 
be an easy miatter to convert the simulation code for the 
controller design into development code for a'^tual 
implementation into hardware. 

For ease of use the programs are menu driven and use 
simple cursor arrows or commanas for pointing at selections. 
Additionally very basic information is needed such as the 
name of functions and equations in a simple format the 
computer can use. This procedure is accomplishea by using the 
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database program which only allows a user to enter 
information in very specific fields of very specific type. 
The type of information required from the designer is to add 
or delete equations, to change initial conditions and to 
provide function names. 

There is an initial learning curve with this program but 
an attempt has been made to keep it as short as possicle. 
There are some peculiarities put in as a requirement to make 
interfacing between the human ana the computer possibl®. Some 
examples include names such as IN?UT:A or IMPUT:E that 
indicate which function the designer is talking about to the 
computer. As an example INPUT :A will associate a function in 
the computer with the function name the designer has selected 
such as INPUT :A = Un . Another very abstract nam= is AlP2r. 
The user is expected to enter either 1 if he wants the 
PRESENT value of the function named under INPUT :A or 2 if he 
wants the EELAYEE value of the function named under INPUT :A. 
If the designer will be using both the present and past 
values then enter 1 to get the present value and the oast 
value will also be available. The strange reouiremer.t fcr 
A1P2D is provided only to allow much faster programs when th^^ 
code is implemented in hardware and the designer only needs 
past or delayed values of a function. Two other codes are lAP 
and lAE for INPUT A PRESENT and INPUT A EELAYEE. These two 
codes are provided to keep the equations needed for the 
functions short. Another variable is ”o” which is usee to 
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define the delayed or previous time sample of the functions 

own OUTPUT. There are similar definitions for up to five 

# 

inputs INPUTiA through INPUT:E for each function. A function 
can be defined in three ways.* means the function is an 
equation that the computer finds a register to hold the 
result in, "AL(x)** means the value is to be placed in one of 
eight analog to digital registers, and '*EA(x)'‘ means the 
value is to be placed in one of eight digital to analog 
registers. An example equation and simple block diagra'T is 
shown in Figure 1. 

Other variables used include ”t" for the discrete sample 
period and "Tl" for the sample period of the plant. If the 
value of T is normalized and Tl has the same value as T then 
the programs will run faster in the realization of the code 
since these values effectively disappear from the equations. 
It should be noted that all functions are done in floating 
point so the simpler the equations the faster the code. The 
example problem at the end of this paper normalizes T but 
still includes it in the equations to show how it could show 
up in the function. Finally the designer is allowed to 
specify end criteria and intervals for calculating the plant 
relative to the computer itself. 

After the designer has entered the equations into the 
database and is satisfied that all function relationships are 
correct then all that is required is to enter the correct 
menu selection that compiles the code. The comnuter program 
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NOTES ; 

Assume by previous definition that INPUT:A = X, INPUT:B = f, 
FUNCNAME = Dz , and the below equation is desired: 

D(k) = X(k) - (2 ♦ X(k-D) + Y(k) - (0.5 * D(k-l)) 

Making substitutions; 

F = Ez = E(k), lAF = X = X(k), lAD = X(k-l), IBP = T = Y(k), 

and 0 = OUTPUT BELAYED = D(k-l) leads to the equation 

F = lAP - (2 * lAD) + IBP - (0.5 O) 

Notice that the above function only has two inputs, X ann Y 



Figure 1 

Simple Difference Equation remonstrating System Variables 
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will now write the necessary source code from the database 
and subsequently invoice the PL/1 compiler which will turn the 
code into a program that will run on any of the previously 
discussed CP/M based systems. The executable file produced is 
called ANSWER.COM and is run on a CP/M based system by typing 
"answer" in response to the operating system input prompt. 
Luring program ANSWER.COM execution the results will be sent 
to a floating point file called LATA.ELT. Also, while the 
program is running the intermediate results are printed on 
the display device for program troubleshooting should the 
results not be as expected. The file EATA.FLT is not 
printable; however, an example utility program called 
PRINT. PLI is provided that changes the floating point file 
into a printable text file. This program is executed by 
typing PRINT FILE where FILE is the name of the designer 
provided text file the designer wants the results placed 
into. With further programming development the floating point 
data can be converted into any format needed by different 
devices such as plotters as shown by the plot results in 
Appendix C. 

The information for understanding the database program 
operation is provided in Chapter V. This section is not 
necessary for understanding at the design level but it does 
explain variables that will be needed in equation 
development. Additionally the example problem will be helpful 
to the design engineer, when the designer is ready to use the 
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to the 



program all that is needed is to enter ”EEA3E MENU 
CP/M operating system prompt. As previously stated the 
designer will need some familiarity with CP/M before the 
program can be of use. 

B. PROGRAMMER INFORMATION 

This section covers the experience necessary for program 
maintenance and further development. These programs would be 
worked on and coded by an individual that has a general 
control system background; however, he must have substantial 
experience with computer programming. The type of computer 
experience needed pertains to operating systems , event 
counters, and automatic programming. The programmer must have 
a good knowledge of both PL/1-S0 and LEASE II. Additionally 
it would be necessary to know microcomputer hardware 
implementations such as analog to digital converters and 
sample and holds. These devices are simulated in the computer 
programs but further useful extensions would involve 
implementing the code into realized hardware. The information 
required for the programmer to understand the programs is 
contained in the remainder of this tnesis. The majority of 
this thesis is the computer code required to implement a 
design environment for a novice programmer through the use of 
automatic programming. 
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III. SYSTEM IIIOSYNCRASIES 



Major problems encountered in tnis thesis had to do with 
programming and efficiency. The major problem is that there 
is no one programming language that contains all the features 
necessary to implement the entire project. EBAS2 II was 
chosen because of the relational database characteristics 
even though it has very limited math or hardware 
capabilities. The requirements for reasonable math is the 
reason PL/l-30 was chosen and though some database systems 
are avaliable for PL/1 they do not have all the features of 
IEAS3 II. A major deficiency of LEASE II is that it will not 
allow the printing of a semicolon. The semicolon is 
considered a command character in LEASE II and therefore 
reserved. This can be considered a major bug for a database 
system, especially one that was designed for business 
applications. It is foolish to oelieve business will not want 
to include a semicolon someplace in one of their 
transactions. Ashton Tate was notified and they replied it 
was not worth trying to fix. This presents a problem since 
the database program was used to write the source code in 
FL/1 and PL/1 requires a semicolon at the end of every 
statement. The problem was worked around by having XBASE II 
write a backslash every place a semicolon was needed in tne 
PL/1 source text. Next a FL/1 preprocessor called REECR'^.PLI 
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was used to change all backslashes to semicolons prior to 
invoking the PL/1 compiler. This takes a large amount of the 
compile time which could have been easily corrected by a 
reasonable response from Ashton Tate. PL/1 also has some 
problems; however, Eigital Research was more than happy to 
try to overcome deficiencies. There are several problems in 
their floating point routines due mainly to conversion r’ules 
which can cause strange effects in the results. Most floating 
point functions such as the exponential function tend to 
convert the arguments to fixed before calculation and this in 
effect truncates the arguments. As an example if one desires 
to raise an exponent to a power and that power is a function 
such as X/Y then both X and Y and the result get truncated. 
If X/Y is a fraction less than 1.0 the effect is to truncate 
the result to 0.0 and tne result of the exponent to a 
fractional power always ends up equal to 1.0. The way around 
this problem is to convert X and Y to CKIHACTERS before doing 
the division. This is only part (f the problem, Even if X and 
Y origionally were floating point and were convertea to 
characters the result of the exponent evaluation is ^’ixed 
decimal. >jow this result must oe converted by anctner 
function called PTC in order to get the result back to tr.e 
system in floating point. An example of this shows up in the 
MONITOR. PLI program for the function callea CP. This does 
cause some limitation on the type of equations that the 
designer can use without knowing some peculiarities of the 



designer can use without knowing sone peculiarities of the 
sys tem . 

For the sake of program efficiency, as far as development 
goes, all variables or functions the designer uses are 
considered floating point. This makes programming easier with 
higher precision; however,, the code takes longer to execute 
even though the execution time was not considered critical in 
this implementation. Using floating point for all functions 
introduces some inaccurate realizations when dealing with 
analog to digital and digital to analog converters. This is 
not a major problem since these devices will have greater 
accuracy in the future or the devices can be modeled mere 
accurately in future designs. 

All results presently produced by this system are placed 
on floppy disks inside the computer, iriting to the disk and 
displaying the results on a terminal are an obvious slowdown 
in system operation. These inefficiencies do not show up in a 
realization of the control system since the controller will 
be self contained computer code communicating only with the 
converters. Additionally the code will be faster since it 
will not be necessary to simulate the plant. 

Presently the realization of the code is not imple'^en ted . 
The software code needed for converters and the actual 
monitor needed for the realized hardware is not developed. 
This is not a major problem and should be relatively straight 
forward for a computer programmer with general hardware 
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experience. The actual control program will he identical and 
will only need to link into the new converter programs. In 
addition the simulation display is not fully implemented. To 
make this project a powerful tool the simulation results 
should he presented in graphical form to the designer similar 
to the plots at the end of this paper. If properly displayed, 
when the designer changes one part of the system all effects 
throughout the system should he dynamically evident. The 
advantages of a completely contained hardware, software and 
display system are obvious. 

The time required to compile the programs is definitely a 
function of the computer hardware. Most newer ZSZ based 
single user systems can compile the code in about fifteen 
minutes. S080 based systems, older ZS0 based systems, and 
multiuser systems that share the central processor will 
require up to thirty minutes to compile. Recall the •'ompile 
time is not critical as mentioned earlier and the realization 
code will run in the order of milliseconds. As a final note 
the programs would run much faster if only one language and 
no preprocessor were needed. 

Another area that requires discussion is the tneory 
involved in the design of the correction factor. levelopmeni 
of a general prediction factor due to delays caused by 
conversions and computer calculations is covered very briefly 
in this paper. The practical results of a predictor that is 
used with the example problem is shown in the graphs at the 
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end of this paper? however, how the results were obtained and 
what they really mean is not fully investigated. Hopefully, 
this will be pursued and discussed in a following paper. The 
predictor is more a control engineers problem and needs to be 
investigated extensively to determine the optimum predictor 
for all possible delays. Further discussion is provided in 
the example problem. 



IV. FUTUR2 EXTENSIONS 



After a careful investigation of the types of problems 
encountered with various hardware and software interfaces, it 
becomes very evident that a sixteen bit or larger 
microprocessor that offers high speed calculation integrated 
into a complete graphic and plotting machine will be a 
definite boost to the realization of future work. 
Additionally a math processor integrated circuit that can 
operate concurrently with the central processor will be of 
invaluable importance. Once you place these hardware items 
along with good converter circuits into one unit it will no 
longer be necessary to simulate in one machine and realize in 
another. Once the simulation code works all that will be 
needed will be to specify executing the realization code on 
the same machine if connectors are available on the machine 
for communicating with the plant. Once the system onerates 
correctly the hardware woula be reduced to a smaller special 
purpose system; however, the power of using a microprocessor 
for the controller cannot be overlooked. If one proviaes the 
capability of a full blown microcomputer in the control loop 
even though all processes are not needed, the system could be 
sent to a distant changing environment and if the control 
equations needed to be changed due to unforseen c i rcums tdnc=s 
it could be done remotely. 2y using the above method it would 
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be possible to reconfigure and realize a new system 
dynamically which cannot be done by a physical one time 
hardware implementation. 

Another requirement is the software system needed for 
this type of development. The problems of integrating two 
completely different software packages is evident in this 
paper. An ideal software package would be a scientific 
relational database that has the capaoility to control the 
computer hardware and provide self modifying code or 
automatic programming for it's own operating system. This 
type of language would make the programming more efficient, 
ease the implementation, and provide a friendlier environment 
for future expansion. This type of software package is 
realizable under the new concurrent languages being developed 
such as ALA (product of the Lepartment of lefence). ALA has 
not been proven to work in microcomputers since it is such an 
extensive language; however, several companies are showing it 
to be feasible in the newer computers such as the 68000 
microprocessor from Motorola. The 66000 is a thirty tw j bit 
processor with a sixteen bit external bus. Motorola is 
presently designing a version of the 66000 to nave a thirty 
two bit external bus. 

It is recommended that future systems contain processors 
such as the 68000 and math processors such as Intel's 3067 or 
similar math processors from other companies such as Tex;as 
Instruments. Additionally high speed high precision analog to 
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dis'ital and digital to analog converters should be used. The 
type of converters used should contain their own multiple 
input sample and hold circuits. This would insure capturing 
all inputs at the same time (versus time multiplexing) making 
the control realizations more practical for real time 
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V. DESCRIPTION OE DATABASE ALGORITHM 



A. GENERAL DISCUSSION (Refer to Appendix A) 

This section covers how the database turns the functions 
into PL/1-80 code. It covers each suoroutiae and there a rp 
hints on how the control designer uses the systerr,. The 
easiest way for the control aesigner to use this thesis is to 
briefly read this material, go through the example pro'^lem. 
and then proceed directly into the program. In general a 
aatabase is a program that contains related information. 
Normally these programs are used in business applications 
where it is desired to keep track of an individuals name, 
address, or account information; however, databases can be 
used to store other information. In this program 
implementation the database contains the fun.ction name, the 
type of function such as an analog function, names of other 
functions that are used by this function as inputs. the 
equations for the functions, the relationship between 
functions, and if the function will be saved for later recall 
and plotting. 

The program has the ability to develop new or recall old 
system models that have been previously aeveloped. Cnee the 
program is entered the model can be changed by aiding or 
deleting transfer functions, change relationships, or change 
the equations. Additionally the models initial conditions can 



be changed. If only the initial conditions are changed the 
program will compile much faster since each function will not 
need to be recompiled. After the system model is developed 
the database proceeds to write PL/1 "include” files that are 
later included into standard PL/1 modules developed for this 
project . 

The features of the LEASE II program are the high level 
relational database language and the ability of the program 
to write text files onto floppy disk. The disadvantages 
include a limited math capability and it does not provide the 
ability to interface with the microcomputer hardware. A 
severe drawback in the program is that it cannot produce a 

semicolon in any of the text files tnat it writes. 

B. LEVEL #1 (MAJOR SUBROUTINES AND PROCEDURES) 

1. MSNy_.CME - initial menu selection 

a. Initializes all variables and the database. 

b. Enters a CASE statement to select a new problem, 

old problem, save the present problem, display present 

function names, change only the initial conditions of the 

present problem, or quit the program. 

c. See Figures 15A and 15B in Appendix A. 

2. CMENy_.CMD - second level menu selection 

a. Develops the control model database. 

b. Enters a CASE statement to edit, add, delete, 
recall deleted functions, or change initial conditions. 
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c 



Calls CONTROL. CML 



d . See Figures 16A and 16E in Appendix A. 

3. CONTROL_.CMr - database control program 

a. Calls procedures BCLMON.CME, FHOC.CML, DATA.C^T. 
and FCRM.CME that develop the FL/1 "include" source code. 

b. Ends the database program and invokes the ?L/1 
programs . 

~ monitor declaration program 

a. Develops DCLKON.TXT for the ['MONITOR .PLI program. 

b. Calls MCNECL.CMr. 

5. ^ONDCL^CMD - develops AA.TXT, AD. TXT, ID. TXT, and 
DA. TXT include functions for '^ONITOR .PLI . 

6- FBQC^CMI - Develops ICLFUN.XXX, INPUT. XXX, FUNC.XXX, 
and OUT. XXX for each function developed by the control 
engineer to be used in the appropriate function AA.PLI, 
AD.PLI, DD.PLI, or DA.PLI. 

P. DATA_.G^E - data declaration proceduf'e 

a. Develops DATAINFO.TXT to be used in PRINT. ?LI 

b. DATAINFO.TXT contains valuable format information 
pertaining to the DATA.FLT file. 

c. Calls DATA1.CMD - additional data format file 

S* - PL/1 command execution program 

a. Develous submit file FORf^.SUP used to control the 
order of PL/1 program compiling. 

b. Issues commands to reformat files containing a 
backslash into files with a sem.icolon. 



c. Calls ECLFUN.CME - function name setup for RtFORi'^ 

d. Calls ESFORM.CMD - general format for PL/1 
compiling. 

9. SHORTiCME - short compile program 

a. Eevelops initial conditions for ECLMON.TXT and 
MONITOR. PLI . 

b. Issues commands to control the compiling of only 
MONITOR. PLI by developing a submit file SFORM.CMI. 

c. Calls SFORM.CME - short format for PL/1 compiling 
similar to FORM. SUB called SFORM.SUE. 

C. LEVEL #2 (PROGRAM DESCRIPTION) 

Since the majority of the code is in higher level 
languages, and is designed in modules, it is already ir 
algorithmic form. It is assumed that the individual needing 
to understand the code will be familiar with references 1 
through 4; therefore, an extensive discussion of the computer 
code will not be provided. The example problem will cover all 
the features of the program. 

The basic structure of the database includes a memory 
file that contains the value of initial variables and a 
database file that contains the functions and their 

relationships to each other. The entire database program is 
involved with manipulating and recalling information from the 
two above mentioned files. 

The initial memory file is shown in Figure 2 and is 
contained in file INIT.MEM. This file is provided so that 
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initial variables called will not produce errors in the 
database program. The main variables of concern are: 

1. T - discrete sample period, 

2. T1 - analog sample period, 

3* PLlljil ~ end criteria for the digital simulation, 

ALIMIT - number of analog samples for each digital 

sample. 

5. E2LAY - delays due to calculation and conversion 
expressed as the number of analog counts after the sample 
period T. 

Unseen system variables include; 

1. ACTION - users next menu selection, 

2* PARAM - holds database program name, 

3. CN - test variable, 

4. FN - test variable, 

5. TEMPOUT - temporary variaole, 

6. ELSMENTS - data array position for a variable. 

7. IND - index counter, 

a. FUN - string of characters that is used to uniquely 
define a function filetyne, 

9* £4141 - number of data array elements saved. 

The basic structure for all functions is shown in Figure 
3 and is contained in file INIT.LBF. The initial values of 
the database are shown at the bottom of the Figure 3. The 
only items the control designer sees are listed under the 
NAME column of Figure 3. The user must enter values for the 



31 



ACTION 


(C) 


0 


PARAM 


(C) 


INIT 


CN 


(C) 


INIT 


FN 


(C ) 


FUNCTION 


T 


(N) 


1 .000 


ELIMIT 


(N) 


10 


ALIMIT 


(N ) 


5 


TEKPOUT 


(C) 


000 .0000 


ELEMEMTS 


(N ) 


0 


IND 


(N) 


0 


FUN 


(c) 


AiCLEFGH 


EELAY 


(N) 


0 


T1 


(N) 


1 .000 


EATAl 


(N) 


0 



Figure 2 

Variables Used for Initialization 
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STRUCTURE FOR FILE: INIT.EBF 

NUMBER OF RECORES: 00001 

BATE OF LAST UPDATE: 00/00/00 



PRIMARY 


USE EAT ABASE 






FLE 


NAME 


TYPE 


WIDTH 


001 


FUNCNAME 


C 


008 


002 


FUNCTYPE 


C 


002 


003 


NR:INPUTS 


c 


001 


004 


INPUT :A 


c 


008 


005 


A1P2L 


c 


001 


006 


INPUT :B 


c 


008 


007 


B1P2L 


c 


001 


00S 


INPUT :C 


c 


008 


009 


C1P2L 


c 


001 


010 


INPUT :E 


c 


008 


011 


D1P2L 


c 


001 


012 


INPUT :E 


c 


008 


013 


E1P2E 


c 


001 


014 


FK 


N 


008 


015 


FEE 


N 


008 


016 


EQUATION 


c 


000 


017 


SAVEEATA 


L 


001 


018 


EATANAME 


C 


005 



EEC 



004 

004 



Below is a list of the values used in the initialization. 
These values correspond directly to the structure above. 

00001 NSW FUNC 
AE 
0 

ONE 

1 

TWO 

1 

TKHEE 

1 

FOUR 

1 

FIVE 

1 

0 .0000 
0 .0000 
AE(1)=0 

.F. 

AD(1) 



Figure 3 
Eat abase Values 
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Ini t ia 1 



variables under NAME of the correct type and size specified. 
The definitions for the variables are given below: 

1. FUNCNAME - function name used to identify this 
function initialized to "NEW ?UNC". 

2. mNCTYPE ~ describes this function's type as; 

a. AA - analog, 

b. AD - analog to digital, 

c. EE - digital, 

d. DA - digital to analog. 

~ This is the number of other functions 
this function uses for inputs. AL(x) and EA(x) registers are 
not counted as functions. 

4. INPUT^A - This is a variable for name of function 
used as an input and is shown at tne bottom of Figure 3 as 
"one". This value is referred to in the equations as lAP or 

lAD for INPUT A PRESENT or INPUT A DELAYED. 

A1P2D - This is a code to tell tne computer the 
designer is talking about the ’*a' input and he wants either a 
"l" for the PRESENT value or "2" for the EELAYSE value. 
Entering a ”l' will make both the present and delayed values 
available but the program will run slower than if only the 
aelayed value is needed. 

6* i^PUT^B - This variable is identical to 4 above 
except this is the "B" input and has an initial name of Ti/O" 
as shown at the bottom of Figure 3. INPUT-.C, INPUTrE, and 
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INPUT:E are the same as above which allows for up to five 
inputs . 

7. B1P2E - This variable is the same as 5 above except 
it is for the ’"b” input function. C1P2E, D1P2D, and E1P2D all 
have the same function except they are used with their 
respective inputs. 

8. FK - This function is the present value of the 
function named under FUNCNAME above. 

S* “ this is the delayed or previous sample period 

value of the function named in FUNCNAME ',ie. FK delayed). 

10. EQUATION - this is the equation for the function 
named in FUNCNAME above. The equation takes three forms! 



F this is a function, 

AE(x) = value for an analog to digital, 

DA(x) = value for a digital to analog. 



11. SAVEEATA - logical value indicating if the user wants 
to save the value of this function for tabulation. 

12. EATANAME - if it is decided to save the function then 
what type function is it! 

F - the function itself, 

AE(x) - one of the eight analog to digital, 

EA(x) - one of the eight digital to anaitg. 

It should be noted that functions are of three types. If 
the computer finds a location for the function or it is not 
hardware dependent then the function is described as ’ ? . If 
the function is one of the other two types "ai(x) or ' In :c)" 
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then they represent physical fixed registers within the 
computer. If the latter two types are used as arguments in an 
equation they are not listed as inputs under MRrINPUTS for 
the function. Since they are physical registers they are 
called by their name (ie. AE(5)) and not by lAP or lAE. 

D. LOGICAL STRUCTaB2 

The logical structure discussed will be a general outline 
of program flow. A complete discussion is covered in the 
example problem. Initially the control designer enters the 
program by responding "IBASE t^ENU" to the CP/f^i operating 
system prompt. This will place the program unaer control of 
MEi^U.CML which initializes the program and enters a CASE 
statement that provides a menu selection shown in Figure 15A. 
At this point the designer usually starts a new problem or 
recalls an old problem that modifications will be made to. 
Notice that in order to save a problem the user must select 
from this menu option. In general the next choice of the 
designer is to enter the modification phase either selecting 
to change only the initial conditions of an old problem or 
selecting to make major changes. If major or new changes are 
chosen the program transfers control to ChiSNU .CtiE where the 
designer is allowed to change the function database by menu 
selections shown in Figure 16A. The choices are controlled by 
the CASS statement in C'^'ENU.CME. ^ost selections are single 
command entries followed by a "RETURN" from the keyboard. 
Some commands will operate without a "RETURN" especially if 
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the user has exceeded a field in one of the database records. 
The fields are marked on the display with a at each end 
of the field. On some terminals the arrow keys may not always 
work for positioning the cursor in a field desired. In this 
case use "CTRL e" (control key pressed at the same time as 
the "e" key) to move the cursor up a field. Use "CTRL X to 
move down a field and "CTRL C" to move down a complete 
reco rd . 

Once the user is sure the database contains all "correct’ 
functions and relationships he should return to the KAIN 
program" if he wants to save the information otherwise any 
changes will be lost during compiling. The program could be 
improved by adding a question in the CONTROL. CME procedure 
asking the user if he wants to save the information before 
compiling. It takes more than fifteen minutes to compile the 
program so a few extra minutes checking for correctness 
before compiling can be well worth the effort. 

Once the option to compile the program is taken, control 
passes to CONTROL. CML. This is the main control procedure 
that causes the automatic programming by calling all 
necessary procedures required to write the PL/ 1 source code 
from the database. A complete discussion of this procedure is 
contained in the example problem. After the PL/1 source .ode 
is written, a control program is written on how to compile 
the PL/1 code. Finally the database transfers control to the 
control program it just wrote, callea FORM. SUB, by exiting 
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the database and submitting FORM. SUB to the CF/^! operating 
system. At this point the PL/l programs are compiled and 
executed . 
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VI. EESCRIPTION CF PL/1 ALGORITHM, 



A. GENERAL DISCUSSION (Refer to Appendix B) 

This section describes the PL/1 programs and how they 
interact to provide simulation of a computer controller for 
an analog plant. There are five basic procedures consisting 
of a monitor that insures all functions are calculated at 
least once during each sample period of the computer, an 
analog proceaure for describing analog calculations. a 
digital to analog procedure and analog to digital procedure 
for describing conversions, ana a aigital procedure for 
describing computer calculations. There is only one monitor 
but there are as many of the other procedures as needed to 
describe the entire system model for each function. As an 
example if there are two analog to digital equations in the 
problem then there would be two analog to digital procedures, 
one for each equation. The example problem will demonstrate 
the function relationships. The basic procedures require 
info rmat ion for each function and that is orovided by the 
control engineer throueh the database program in the for"^' of 
include files. Once the program is ready and compiled 
automatically by the database program a result file called 
ANSWER.COM is produced. The problem is run by typing "anJWER' 
to the CP/M operating system prompt. As the program is 
running intermediate results are displayed on the display 
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device for troubleshooting should the results not be as 
expected. After the program runs it produces a file DATA.FLT 
containing floating point data. This file can be translated 
into a text file by typing PRINT FILS where FILE is a name of 
a file provided by the user that he wants the results placed 
into . 

PL/l was chosen as the operating system language because 
of several features. First it is a common high level language 
that is transportable to many of the present microcomputer's. 
The intent of this project was lo 'teep the pr'ogram s as 
machine independent as possible and using FL/1 source code is 
a good way to do this. PL/l also has very good math 
capabilities. Since the type of problems being solved are 
engineering in nature it was necessary to have a high level 
math language. Additionally, the ability for PL/l to 
communicate with the computer hardware is a definite asset if 
the programs are to be realized into hardware. Several other 
languages were looked at such as PASCAL but PASCAL does not 
have the ability to talk to the hardware. The 'C language 
was also looked at but it lacks the math capability. As 
stated in Chapter IV.; hopefully, ADA will cverccme all the 
difficulties by combining all the needed features into one 
language . 

PL/1-S0 has some problems. There is a definite nroolem 
with conversion factors as was discussed in Chapter III. The 
main problems have to do with conversions between floating 
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point and other forms. The original approach was to use only 
floating point values and therefore not need conversions but 
it was later discovered that several functions do conversions 
regardless of starting and ending types. This is a problem 
and cannot be overlooked by the designer if he has problems 
with results. If equations requiring complex functions are 
used the designer will need to become familiar with PL/1. 

L. LSVZL #1 (r^.UOH SUIR0UTIN23 AND PLOCSIjHIS) 

1. i^ONITOR.PLI - main control loon urogram 

a. Causes iteration of ail functions. 

b. Insures all functions are callec in proper order. 

c. Saves data for later display and analysis. 

2. AA_.?LI^ - basic format for all analog functions 

a. Calls all functions used as inputs. 

b. Calculates and returns it's own present or oast 
values to all functions that call it. 

c. Insures synchronization of calculations with the 
system clocks through event counts. 

“ basic format for all analog to dic-ital 
functions and performs identical functions to those shown in 
2a, 2b, and 2c above. 

4:. Tli?LI - basic format for all digital functions and 
performs identical functions to those shown in 2 above. 

" basic format for all digital to analog 
functions and performs identical functions to those above. 



6. ~ corrects all improoerly formed PL/1 
"include" files by changing the backslash produced by the 
BDASE II program into semicolons. 

7. PRINT^PLI - utility program showing how to convert the 
floating point file EATA.PLT produced by the main program 
into a different printable text file for display on the 
system display. 

C. LEVEL #2 (PROGRAM EESCRIPTION) 

Since the majority of the code is in higher level 
languages, and is designed in modules, it is already in 
algorithmic form. It is assumed that the individual needing 
to understand the code will be familiar with references 1 
through 4; therefore, an extensive discussion of the computer 
code will not be provided. The example problem will cover all 
the features of the program. 

The basic structure of the program uses a monitor that 
sequentially calls all functions within the computer once 
during each digital clock. The system clock is used as an 
event counter to provide synchronization for each function. 
The functions are developed by the control engineer through 
the database program and each functional relationship 
indicating which functions call each other must be nrooerly 
provided. This is taken care of when the designer states 
which functions are inputs to each other. A function can have 
its output connected to any number of functions but a 
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function can only have five inputs not counting analog to 
digital, digital to analog, and it's own delayed output. The 
function definition allows for very complex multiple 
input/output relationships within any single equation. A 
method for correcting for computation and other delays is 
provided in the correction factor CF in the monitor. This is 
a factor provided for the example problem; however, it is 
felt that this equation will work for any control system of 
the type requiring a steady state error signal of zero. The 
CF is provided hut tan be eliminated simply by not 
referencing it in any of the equations. The programs provided 
are for simulation only but should be easily realizefi into 

control hardware since the design equations will not change. 

# 

L. LOGICAL STRUCTURE 

The program contains a numbe r of procedures that are 
basic to all problems. These were discussed above. 
Additionally each problem has unique features that are 
provided by the control designer through the database 
programs. These features are coded into PL/1 "include” files 
and read in at compile time. Once the programs are coded all 
information as to which functions are connectec together will 
be provided. The monitor is a DO a^HILS LOOP that calls all 
analog to digital and digital functions once each count of 
the clock. At each count of the digital clock ( DCLOCK ; any 
called function will perform it's own equation evaluation if 
it has not already been calculated during this cIock and 
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return the correct value to the invoking function. When the 
monitor calls a function it does not use the value returned. 
Additionally there is a DO WHILE LOOP inside the main monitor 
loop used to call the analog functions the number of times 
specified by the designer to occur between each digital 
sample period. This is the analog clock ( ACLOCK ). Also 
during the analog clock when the ACLOCK is equal to the DELAY 
specified by the designer the digital to analog function is 
invoked. When a function is invoked, the first thing that it 
does is check Co see if it has any inputs. If there are any 
inputs then the input functions are called to insure they are 
the correct value before the equation is calculated. Sven 
though the monitor is a sequential process invoking each 
function, once a function is called it will invoke all of the 



input functions 


needed for equation 


calculation. 


Some 


functions could 


be called 


many times 


luring each 


clock 


depending on how 


many other 


functions use 


it's output . 


Once 


each ECLOCK the 


monitor 


outputs values that have 


been 



requested by the designer to save into a file. 

The first procedures called by the monitor are all 
digital to analog functions located in a file written by the 
database program called "ad.txt". The functions take the form 
iummy_ptr = ?UNCTI0N_NAMS where EUNCT I0 N_i'jAMS is the name of 
the function and the file EUiMCT I0N_NAME contains the 
necessary information to represent the function. The example 
problem has a more detailed discussion. The first thing the 



44 



digital to analog functions do is call the analog functions 
that 'are used as inputs. If these analog functions have other 
analog inputs they call their inputs before calculating their 
output results. This process will continue until the analog 
functions have no inputs or their inputs come from a digital 
to analog register that will contain a value for a correct 
point in time. Similar processes occur for files "dd.txt", 
"da.txt", and ’’aa.tit" representing digital, digital to 
analog, and analog functions respectively. The crder of 
calling is controlled initially by the monitor and is 
extremely important for all processes to obtain the correct 
results. Direct coding of the processes would be much more 
efficient but it must be remembered that tnese processes must 
work no matter how- the designer initially defined the 
problem. This method of problem solution has merits even for 
the experienced programmer. Even though the code is longer 
than direct programming this type of problem will typically 
require one or two hundred files introducing the possibility 
for syntax and other programming errors. 

The four main functions aa.pli, ad.pli, dd.pli, and 
da.pli all have the same general structure. Easically the 
function is called by the monitor cr another function, the 
function calls any inputs, and then returns the correct value 
of the equation at the time the equation was invoked. This 
arrangement places a process on the stack through the FL /I 
calling routine until all inputs are available eliminating 



the need for developing a special process nianager. Using this 
approach makes the programs easily adapted to concurrent 
processors and languages such as AEA. The conversion routines 
ad.pli and da.pli do not have a provision for using their 
delayed values since they would not be available in real 
hardware. This forces a true representation of the hardware 
and makes the designer provide intermediate holding functions 
if delayed values are needed. Additionally the actual code 
for hardware will be easier to implement at a later date. The 
other functions dd.pli and aa.pli have provision to return a 
previously calculated value if only the delayed value is 
desired. This allows the code to run much faster and provides 
more powerful equations within a single function. All four 



processes 


allow for 


up 


to five function 


inputs and 


the 


procedure 


knows how 


to 


call those inputs 


by using 


the 


variable 


"nr_inputs " 


to 


determine to which 


label in 


the 



function to jump, basically all functions the designer 
implements are written into one of the above four procedures. 
As an example if the designer specifies a digital function 
called DIGIT, the database will write the include files and 
then compile the dd.pli procedure into a file DE.REL. Next 
the compiled file ED.RSL will be copiea into DIGIT. REL wnich 
will now contain all the necessary information for the 
function DIGIT. Additionally the necessary information such 
as the name DIGIT will be included into the monitor so that 
it knows to call the function once each clock period. 
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Before any of the PL/1 programs can be compiled, all the 
text files written by EBASE II that contain a backslash must 
be converted to semicolons. This is accomplished by the 
procedure REFORM. PLI. This is a simple program that reads in 
a file, translates the character, and outputs to a new file. 
The files reading from and too are controlled by the submit 
files FORM. SUB or SFORM.3UB written by the EBASE II program. 

The programs compile into a program called NS'i/EH .COM 
under the control of the database program. ANSWER. COM is the 
file needed to run the problem. The designer only needs to 
place this file on a CP/M based system and type AIvSWSr! in 
response to the operating system prompt. This will run the 

w 

program and sand the intermediate results to the display and 
send the floating point data into a file called DATA.FLT. 

Once the program has been run it may be desired to nlace 
the floating point data into a different form for display and 
analysis. This project can be greatly improved oy providing 
the results on a graphical display device. Both tacular and 
graphical results are supplied at the end of this paper and 
the reader can determine for himself which is the best 
m.ethod. One utility program provided is called PnlNT PLI. 
This program reads in the floating point data by using the 
information contained in the file DATAINFO.TXT as to how the 
data was stored and then sends the information into a text 
file specified by the user. This text file can then be 
printed if desired. 
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VII. EXAMPLE 



A. INTRODUCTION 

This part of the paper will make reference to all 
aspects of the project and provides the easiest method for 
understanding the entire paper. It should be noted that the 
descriptions as to how the data is entered and displayed is 
somewhat misrepresented since the entire project was designed 
for an interactive environment. The actual data entry at 
runtime is much more formal with data fields that allow only 
specific entry in specific areas of very specific t'ype . This 
is not a constraint but an aid. As an example if the user is 



to 


replay 


if he wants 


to save the 


function 


resul ts 


he 


can 


only 


answer 


with ”T , ? 


, Y , or N" 


for "true 


, false. 


yes, or 


ft 

no . 


The 


answers can 


he in upper 


or lower 


case and 




an 



improper input is made such as "7” the computer will BEEP and 
wait for a correct response. Additionally some information is 
supplied in reverse video or highlighted. In the event that a 
color system is used such as the INTERCOLOR the information 
will be in different colors. The only way to fully appreciate 
the aids for the designer is to run the program. In order to 
provide a realistic test for this program the sample problem 
will include a plant that is at the limit of stability and is 
initially designed to be controlled by a bang bang 
controller. The process will include a control loop that 



feeds back on itself which provides a definite headache to 
computer programmers because of the synchronization problems. 

b. bKOkLKM SKTUP 

2 

Assume we have a continuous 1/3 plant and it is desired 
to produce a minimum time response ripple free controller 
L(k) for a unit step input. Referring to Figure 4A below, 
L(k) is the computer transfer function, T is the sample rate 
of the input, E is the delay from input sample at time T to 
output sample including all computation and conversion 
delays, and the clamp is included since the computer output 
is held fixed between samples. Assuming that the clamp has a 



T T + D 




Figure 4A 

Computer and Plant Model 




very short conversion delay compared to all other delays in 
the digital system, the output of the clamp is for all 
practical purposes the same as the input to the clamp. This 
has been demonstrated by many authors [Eef. 5]. The next 
initial assumption will he to ignore any delays due to 
conversions and calculations in the computer Ue set E = 0). 
This will allow for simple design equations yet the gross 
error in reasoning will be easily corrected for later. ' It 
should be noted that most control systems ignore the delays 
due to computer computation since the delays are usually 
insignificant compared to the sample rate. ) 

By letting X2(t) = the plant .output , Xl(t) = d/dt[X2U^]. 
and letting the input to the plant be E2(t) by ignoring the 
effects of the clamp and computation delay, a set of discrete 
difference equations can be developed describing the plant at 
the sample instances T. These equations have been developed 
in the past by many authors [Ref. 6] and are shown belrv in 
discrete form: 

Xl(k+1) = Xl(k) + T * E2(k) (1) 
X2(k+1) = T - Xl(k) + X2(k) - [(T T) / 2] * 12 (k; ('2; 

where; 

d/dt[Xl(t)] evaluates to Xl(k+1) 
d/dt[X2(t)l evaluates to X2(k+1) 

Now the next task is to control the plant described by 
the equations in 1 and 2 above according to the design 
requirements for a minimum time response ripple free 
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controller. Since this is a second order system it will take 
two time periods to obtain control. Assume the initial state 
vector equals zero ( ie . X1(0) = X2(0) = 0). Since we want the 
second time state with Xl(2) = 0 or no ripple and X2{2) = 1, 
the input unit step, this leads to the set of equations below 



after some manipulation: 

Xl(2) = T S2(0) + T * S2(l) = 0 (3^ 

X2(2) = [3’MT * T)/2]*E2(0) + [(T * T)/2]*E2(1) = 1 (4) 

Solving equations 3 and 4 above simultaneously leads to* 
E2(0) = 1/(T - T) (5) 

S2(l) = - 1/(T * T) (6) 



By definition of the Z transform, equations 5 and S and 

9 

the fact that the desired output requires E2(k) = 0 for k > 1 
leads to: 

-1 

E2(z) = 1/(T - T) - [1/(T - T)]z ^7) 

Now to find the Z transform of Sl(k) we start with the 
problem definition of zero initial conditions and a unit step 
input . 

At k = 0 

S1(0) = U(0) - X2(0) 

= 1 - 0=1 
A t k = 1 

Sl(l) = U(l) - X2(l) 

= 1 - [(T * T)/2] - S2(0) 

=1 - [(T - T)/2] - 1/(T - T) 

= 0.5 



Cl 



1 



0.5 



At k = 2 

21(2) = U(2) - X2(2) 

= 1 - 1 =0 
For K > 1 EKk) = 0 

The above set of equations leads to the Z transform: 

-1 

EKz) = 1 + 0.5z (8^ 

Since r(z) = 22 ( z ) / El(z), the algorithm for the 
computer, we can now take the ratio of equations 7 and 8 and 
cross multiply then transform to find the difference 
equations needed. 

-1 -1 
[l + 0.5z ]*E2(z) = [l/(T - T)] « [l - z ]*El(z) (5) 

E2(k) + 0.5?E2(k-l) = [1/(T * T)]-[Sl(k) - El(k-l)] (10) 

E2(k) = - 0.5-'!‘E2(k-l) + [l/(T * T ) ] * [El ( k ) -El ( k-1 ) ] 'll) 

Equation 11 is the one used for the computer simulation 

and shown in Figure 4E along with the rest of the model. 

Since it is not allowed to use subscripted or arrays for 

variable names in this program implementation the following 

changes have been made in Figure 4B for variables in 

equations 1, 2, and 11: 

SK = Sl(k) 

Sl(k-l) referenced as the delayed value of EK or lAL in 
the function E2K . 

E2K = E2(k) 

E2(k-l) referenced as the delayed value of E2K or "O in 
the function E2K . 

XIN = Xl(k) or Xl(n) if n different from k 




figure 4B 

Test Program Block Eiagram 
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Xl(k-l) referenced as the delayed value of XIN or "o" in 
the function XIN and lAI in function X2N1. 

X2N1 = X2(k) or X2{n) if n different from k 
X2(k-1) referenced as the delayed value of X2{k) or "O 
in the function X2N1 and lAI in function AD2. 

UN = U(t) 

T1 = [T * i] / j 

In the initial design the sample rate of the digital 
system T and the sample rate of the analog system at discrete 



points 


T1 


were considered 


the same( ie. i 


= j). The 


only 


reason 


for 


providing the 


possibility of 


calculating 


the 


analog 


values at different 


points is to see 


what happens 


to 



the plant between sample points of the digital system at T. 
In most problems T and ”i" will be normalized to unity and 
"j” will represent the number of sample points calculated by 
the plant for each unit sample period T. 

Now to account for the combined delays in the discrete 
system due to analog to digital, digital to analog, and 
computation delays. This function is called CF or the 
Correction Factor. The GF in most systems is a major 
consideration in initial equation development and usually 
resembles Kalman Filters or predictors. The approach taken 
for this thesis was that the average designer would not know 
or could not estimate these delays due to system variables. 
This means the simulations would not work as expected and 
needless to say neither would the realizations of the control 






algorithms in the real hardware. The attempt was to force the 

designer to consider the delays in his design and possibly 

provide a good predictor initially for his use. It was 

assumed that since this type control system requires an error 

signal that is zero when in control this would indicate that 

the signal is zero mean. Since the instantaneous input is 

random it was assumed the error signal was gaussian and 

provide a standard exponential predictor to the error signal 

-I 

to correct for delays. The predictor chosen was CF = e 
where e is the natural log and E is the combined delays. The 
CF is not a complete guess and further work will be presented 
at a later date. The results of the CF to this particular 
example can be seen in the plots at the end of this paper. 
Further research needs to be done to see if this correction 
works for truncation, roundoff, or even noise in the input 
signal. The point is that some type of correction needs to be 
added and this one works very well in this case. If the CF is 
not wanted the designer only needs to eliminate it from any 
of his equations . 

C. AUTOMATIC PROGRAM (Refer to Appendix A and E) 

The designer will initially enter the program by entering 
' LEASE menu" and a RETURN to the CF/M operating system 
prompt. This assumes the disk containing the ESASE programs 
is in drive "a and the disk containing the PL/1 programs is 
in drive ’ E" . The system can be modified to run on a single 
aisk as was done in the original system but the disk must 



contain at least one megabyte and should have directory space 
for at least one hundred fifty files. Ey using the program 
setup requirements above the database is entered by menu 
selection "l" of Figure 15A to start a new program. This will 
initialize the system and enter the menu selection of Figure 
16A. Next if menu selection "s" is entered the initial 
conditions of Figure 2 can be changed to the ones shown in 
Figure 5A. This is done by the display menu CINIT of Figure 
17A. Note that some parameters of Figure 5A are effected by 
returning to the main menu and saving the program under the 
name STRT, others are developed as part of the program 
process such as ACTION which changes each time a menu 
selection is made, and others such as EATAl are developed by 
the program itself counting the number of items that have 
been requested to save. As can be seen from Figure 5A , 
ALIMIT was entered as 13 which means there will be 10 analog 
calculations for each digital calculation. Since the initial 
design used the same sample time for both the digital and 
analog then T1 must be adjusted to 0.1 as discussed 
previously since T is normalized to 1. After the initial 
conditions ALIMIT, DLI^^IT, DELAY, T, and T1 are entered, the 
program returns back to CMENU of Figure 16A. 

Next assume menu selection " 2 " is selected and the 
functions are added. The structure will remain basically the 
same with only slight changes shown in Figure 5?. The 
database for the problem will look something like that shown 
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ACTION 


(C) 


3 


PARAM 


(C) 


STRT 


CN 


(C) 




FN 


(C ) 


FUNCTION 


ELEMENTS 


(N) 


0 


ELIMIT 


(N) 


10 


ALIMIT 


(N) 


10 


T1 


(N) 


0.100 


FUN 


(C) 


ABCEEFGHI JKLMNOPQRSTUVWXYZ 


INL 


(N) 


7 


T 


(N) 


1 .000 


EELAY 


(N) 


4 


EATAl 


(N) 


9 






Figure 5A 



Variables Used in STRT Initialization 



structure for FILE: 


3TR 


T .EBF 


NUMBER 


OF recores : 


00008 


eate of 


LAST UPDATE 


: 11/11/82 


PRIMARY 


USE DATABAS 


E 




FLD 


NAME 


TYPE 


WIDTH 


001 


FUNCNAME 


C 


008 


002 


FUNCTYPE 


C 


002 


003 


NR '.INPUTS 


c 


001 


004 


INPUT ;A 


c 


008 


005 


A1P2E 


c 


001 


006 


INPUT:B 


c 


008 


007 


E1P2E 


c 


001 


008 


INPUT:C 


c 


008 


009 


C1P2E 


c 


001 


010 


INPUT :E 


c 


008 


011 


D1P2D 


c 


001 


012 


INPUT :E 


c 


008 


013 


E1P2E 


c 


001 


014 


FK 


N 


008 


015 


FEE 


N 


008 


016 


EQUATION 


C 


050 


017 


SAVEEATA 


L 


001 


015 


EATANAME 


C 


005 




Figur 


e 5B 






Struc ture 


for 


STRT 
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in Figure 6A - 61 beloi«. Keep in mind that the display looks 
quite different due to color, highlighting, or reverse video. 
Also, only one record at a time will he on the screen. 

The remainder of the menu programs will not he covered 
since they are self explanatory. It should be noted that in 
order to save the problem the designer is working on, he must 
first return to the main menu and then save the parameters. 
The program will ask for the name of the file the user wants 
to save the program in. The example problem was saved in file 
3TRT . If there are questions about the menus look at 
procedures I^ENU.Ct^T and CMENU.CME in Appendix A. After all 
functions are entered CORRECTLY, and the program is saved if 
desired, then enter ” 0 " for menu selection in Figure 16A. 
This will cause the program to enter the automatic 
programming phase. (Take a break because it will take about 
fifteen minutes.) If the program fails to compile it will 
probably be due to improper naming of functions. Function 
names cannot be arrays. If the results are not as expected, 
even though the program compiles, then the equations are 
probably in error or there is a strange conversion in PL/1. 
If this happens then unfortunately references 2 through 4 
must be studied. 

Once compiling begins control passes to CONTROL. CME. The 
first function called is ECLMCN .CME which writes the monitor 
declarations ICLf^ON.TXT for MONITOR. PLI and files A&.TTT, 
AI.TTT, EE.TTT, and EA.TTT containing the procedu’^e calls for 
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RECORD 00001 
FUNCN.\ME 
FUNCTYPE 
NR .‘INPUTS 
INPUT: A 
A1P2E 
INPUT:E 
B1P2D 
INPUT:C 
ClP2r 
INPUT :D 
E1P2D 
INPUT:E 
E1P2D 
FK 
FEE 

EQUATION 

SAVEDATA 

EATANAME 


UN 

AA: 

0 : 

ONE : 

1 : 

TWO : 

1: 

THREE : 

1 : 

FOUR : 

1: 

FIVE : 

1 : 

1 .0000 : 
0 .0000 : 
F=1 .0000 
T : 

F : 


RECORD 00002 
FUNCNAME 
FUNCTYPE 
NR: INPUTS 
INPUT: A 
A1P2D 
INPUT:B 
E1P2E 
INPUT:C 
C1P2D 
I N PUT : E 
D1P2E 
INPUT :S 
E1P2D 
FK 
FKD 

EQUATION 

SAVEDATA 

DATANAME 


EX : 

ED: 

0: 

ONE : 

1 : 

TWO : 

1: 

THREE : 

1: 

FOUR : 

1 : 

FIVE : 

1 : 

0 .0000 : 

0 .0000 : 

F=(AE(l)-Ar^2M*GF 

T : 

F : 




Figure 6A 



Database For Example Problem 

t=.C 



RSCORI 00003 
FUNCNAME 
FUNCTYPE 
NR -.INPUTS 
INPUTtA 
\1P2L 
INPUT.-B 
B1P2D 
INPUT:C 
C1P2E 
iNPUTtD 
E1P2E 
INPUT:E 
E1P2E 
FK 
FO 

EQUATION 

SAVEEATA 

DATANAME 


E2K : 

EE: 

1 : 

SK : 

1 : 

TWO : 

1: 

THREE : 

1 : 

FOUR 

1 : 

FIVE : 

1 : 

0 .0000 : 

0 .0000 : 

F = (IAP-IAE)/(T’:‘T)-.5’^0 

T : 

F : 


RECCRE 00004 
FUNCNAME 
FUNCTYPE 
NR :INPUTS 
INPUTtA 
A1P2E 
INPUT :B 
B1P2E 
INPUT :C 
C1P2E 
INPUT :E 
D1P2E 
INPUT :S 
E1P2E 
FK 
FKE 

EQUATION 

SAVEEATA 

EATANAME 


XIN : 

AA : 

0 : 

ONE : 

1: 

TWO : 

1 : 

THREE : 

1 : 

FOUR : 

1 ; 

FIVE : 

1 : 

0 .0000 : 

3 .0000 : 
F=EA(l)''‘Tl+0 
T : 

F : 




Figure 6B 



Eatabase For Example Problem 
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RECORD 00005 
FUNCNAME 
FUNCTTPE 
NR:INPUTS 
INPUT; A 
A1P2D 
INPUTtB 
B1P2E 
INPUTtC 
C1P2I 
INPUT :E 
D1P2D 
INPUT :E 
E1P2E 
FK 
FKL 

EQUATION 

3AVEDATA 

DATANAME 



X2N1 

AA; 

1 : 

XIN : 

2 ; 

TWO : 

1 : 

THREE : 

1 : 

FOUR ; 

1 : 

FIVE ; 

1 : 

0 . 0000 ; 

0 .0000 : 

F=EA(1 )*T1*T1/2+IAB*T1-^0 
T : 

F 



RECORD 00006 




FUNCNAME 


A El 


FUNCTYPE 


AE: 


NRiINPUTS 


1; 


INPUT:^ 


UN : 


A1P2D 


1 : 


INPUTtB 


TWO 


B1P2B 


1 : 


INPUTtC 


THREE 


C1P2D 


1 : 


INPUT:! 


FOUR : 


D1P2D 


1 ; 


INPUT :E 


FIVE : 


S1P2E 


1 : 


FK 


0 .0000 : 


FKD 


0 .0000 : 


EQUATION 


AE(1)=IAP 


SAVEEATA 


(V • 


DATANAME 


AD(1) : 



Figure 6C 

Database For Example Problem 
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RECORE 00007 


FUNCNAME 


AD2 


FUNCTYPE 


AE; 


NRrINPaTS 


1: 


INPUT: A 


X2N1 


A1P2E 


2: 


INPUT ;B 


TWO 


E1P2E 


1: 


INPUT:C 


THREE 


C1P2E 


l: 


INPUT :E 


FOUR 


E1P2E 


1 : 


INPUT :S 


FIVE 


E1P2E 


1: 


FX 


0 .0000 


FKD 


0 ,0000 


EQUATION 


AE(2)=IA 


SAVEEATA 


T: 


DATANAME 


AD(2) : 


RECORE 00009 


FUNCNAME 


EAl 


FUNCTYPE 


EA: 


NR: INPUTS 


1 ; 


INPUT :A 


E2X 


A1P2E 


1 : 


INPUT:B 


TWO 


B1P2D 


1 : 


INPUT:C 


THREE 


C1P2E 


1 : 


INPUT:E 


FOUR 


E1P2E 


1 : 


INPUT :S 


FIVE 


S1P2D 


1 : 


FK 


0 .0000 


FKE 


0 .0000 


EQUATION 


DA(1)=IA 


SAVEEATA 


O’ . 

1 . • 


EATANAME 


EA( 1 ) : 



Figure 6D 

Eatabase For F.xample I 
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the monitor. See Figure 7. Next PROC.CMI is invoked for each 
of the four types of functions. This procedure writes the 
declarations for each function called CCLFUN.XXX where XXX 
stands for which function number it is (ie. A, B, C etc. and 
the last two XX indicate the function type (ie. kk, AI, IE, 
EA). PROC.CME also writes the function input files INPUT. XXX, 
function equation files FUNG. XXX, and the output files 
OUT. XXX that are used to indicate which element in the data 
array each function represents. See Figures 8A through 8H. 
The next function invoked is EATA.CMD. This procedure writes 
the file EATAINFC.TXT that tells how many functions are 
stored in the floating point file EATA.FLT and the location 
of each element. See Figure 9. Next the program passes 
control to FORM.CME. This procedure writes the program that 
tells the computer what to do with all the files that have 
just been written. Basically it produces a CP/M SUBMIT file 
that tell the system to reform or correct all files that have 
a backslash into files containing semicolons, then it tells 
the system to compile each function and place it under it's 
own unique function name. After compiling each function it 
links them all together into one program called ANSV;ER.C0M. 
See Figure 10. Finally after the database program writes the 
SUBMIT program it passes control to the program it just wrote 
and creates the simulation program ANSa’ER.COM as just 
ment ioned . 
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/« rCLMON.TXT */ 



UN ext entry ( fixed (7 ) ) returns ( ptr ) , 

XIN ext entry ( fixed (7 ) ) returns (ptr ) , 

X2N1 ext entry ( fixed (7) ) r eturns ( ptr ) , 

All ext entry( f ixed (7) ) returns ( ptr ) , 

AI2 ext entry ( fixed (7 ) ) returns ( ptr ) , 

DAI ext entry( f ixed (7) ) returns ( ptr ) , 

EX ext entry (fixed (7 ) ) returns (ptr) , 

E2K ext entry ( fixed (7 ) ) returns ( ptr ) , 

T float static ext init( 1.0000), 

Tl float static ext init( 0.1000), 

DLIMIT fixed static init( 10), 

ALI^’IT float static init( 10), 
delay float static init( 4), 

1 dat ext, 

2 data( 3) dec(15,6) 



/- AA.TXT */ 



durnmy_ptr 


=z 


UN 


( code ) \ 


dummy_ptr 


= 


XIN 


( code )\ 


dummy_pt r 


rr 


X2N1 


( code )\ 


AD. TXT */ 


dummy_pt r 


= 


ADI 


( code ) \ 


dummy_pt r 




AD2 


( code )\ 


DD.TXT */ 


dummy_ptr 


= 


EX 


(code )\ 


dummy_pt r 




E2X 


( code )\ 


DA. TXT -/ 


durnmy_pt r 


= 


DAI 


( c ode ) \ 



Figure 7 

r^onitor Declarations and Functions 
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/* DCLFUN.AAA -/ 



UN :proc(code) returns ( ptr)\ 

del 

nr_iaputs fiied(7) static init(0), 

ONE ext entry ( fixed (7 ) ) returns (ptr ) , 

TWO ext entry( f ixed( 7) ) returns ( ptr ) , 

THREE ext ent ry ( f ixed ( 7 ) ) returns (ptr ) , 

FOUR ext entry( f ixed (7) ) returns (ptr ) , 

FIVE ext ent r y ( f ixed ( 7 ) ) returns (ptr ) , 

cod(5) fixed(7) static init ( 1 ,1 ,1 ,1 ,1 ) , 

F float static init( 1.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data( 8) dec(15,6)\ 

/« ECLINP.TXT - THIS IS PART OF ALL FUNCTIONS ANE WAS 

FIXED FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. -/ 

del 

1 lA oasedUa ptr), 

2 lAP floatT 
2 lAE float, 

1 13 based(id ptr), 

2 ISP floatT 
2 ISE float, 

1 I C based ( ic_ptr ) , 

2 ICP float, 

2 ICE float, 

1 ID based ( id_pt r ) , 

2 lEP float, 

2 lEE float, 

1 IS based ( ie_pt r ) , 

2 ISP float, 

2 lEE float; 

/* INPUT. AAA */ 
input(l) = ONE \ 

input U) = TWO \ 
input(3) = THREE \ 
input(4) = FOUR \ 
input(5) = FIVE \ 

/* FUNC.AAA */ 

F=1.0000 

/-' OUT. AAA -/ 

da ta ( 1 ) = f t c ( F ) \ 

put skip list('data( l) = ',F )\ 



Figure 8A 

Function Declarations and Equations 



/* ECLFUN.BAA */ 



XIN :proc(code) returns (ptr )\ 

del 

nr_inputs fixed(7) static init(0), 

ONE ext ent ry ( f ixed (7 ) ) returns (ptr ) , 

TWO ext entry( f ixed(7) ) returns (ptr ) , 

THREE ext ent ry ( f ixed ( 7 ) ) returns ( ptr ) , 

FOUR ext entry(f ixed(7) ) returns ( ptr ) , 

FIVE ext entry( f ixed(7) ) returns ( ptr ) , 
cod(5) fixed(7) static ini t (1 , 1 , 1 , 1 , 1 ) , 

F float static init( 0.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext , 

2 data( 8) dec( 15,6)\ 

/* CCLINP.TXT - THIS IS PART OF ALL FUNCTIONS AN! WAS 

FIXEE FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. */ 

del 

1 lA oased( ia_ptr ) , 

2 lAP float, 

2 lAE float, 

1 IB based( ib ptr), 

2 IBP floatT 
2 IBE float, 

1 I C based ( ic ptr), 

2 ICP floatT 
2 ICE float , 

1 IE based( id_ptr ) , 

2 lEP float , 

2 lEE float, 

1 IE based( ie_ptr ) , 

2 ISP float , 

2 lEE float; 



/* INPUT. BAA «/ 




input(l) = ONE 


\ 


input(2) = TWO 


\ 


input (3) = THREE 


\ 


input(4) = FOUR 


\ 


input(5) = FIVE 


\ 


/* FUNC.BAA -/ 




F=EA (l)-Tl+0 




/* OUT. BAA -/ 





data( 2) = ftc(F )\ 

put skip list('data( 2) = ",F )\ 

? i gu r e SI 

Function Eeclarations and Equations 



/* ECLFUN.CAA -■!'/ 

X2N1 :proc(code) returns (ptr )\ 
del 

nr_inputs fixed(7) static init(l), 

XIN ext entry( f ixed(7) ) returns (ptr ) , 

TWO ext ent ry ( f iied (7 ) ) returns (ptr ) , 

THREE ext entry( f ixed(7 ) ) returns (ptr ) , 

FOUR ext entry( f ixed( 7) ) returns (ptr ) , 

FIVE ext entry( f ixed (7 ) ) returns ( ptr ) , 
cod(5) fixed(7) static ini t ( 2 , 1 , 1 , 1 , 1 ) , 

F float static init( 0.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data ( 3 ) dec ( 15 ,6) \ 

/* ECLINP.TXT - THIS IS PART OF ALL FUNCTIONS AND WAS 

FIXED FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS . -/ 

del 

1 lA based (ia_ptr ) , 

2 lAP float , 

2 lAD float, 

1 IB based ( ib_pt r ) , 

2 IBP floatT 
2 IBI float, 

1 IC based ( ic_ptr ) , 

2 I CP float, 

2 ICL float, 

1 ID based (id ptr), 

2 lEP floatT 
2 IDD float, 

1 IE ba sed ( ie ptr), 

2 lEP floatT 
2 lEE float; 

/- INPUT.CAA -/ 
input(l) = XIN \ 
input (2) = TWO \ 

input (3) = THREE \ 
input(4:) = FOUR \ 
input(5) = FIVE \ 

/* FUNC.CAA */ 

F=DA(l)-Tl-Tl/2+IAE-Tl+0 

/- OUT.CAA -/ 

data( 3)=ftc(F )\ 

put skip list('data( 3) = ',F )\ 

Figure SC 

Function Declarations and Equations 
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/- DCL?UN.AAD 



All :proc(code) re turns ( pt r ) \ 

del 

nr_inputs fixed(7) static init(l), 

UN ext entry( f ixed(7) ) returns (ptr) , 

TWO ext ent ry( f ixed ( 7 ) ) returns(ptr ) , 

THREE ext ent ry ( f ixed (7 ) ) returns (ptr ) , 

EOUR ext ent ry( f ixed (7 ) ) r eturns ( ptr ) , 

FIVE ext entry (fixed(7 ) ) returns ( p tr ) , 

cod(5) fixed(7) static ini t U , 1 , 1 , 1 ,1 ) , 

F float static init( 0.2000), 

0 float static init( 0.0000)\ 
del 

1 dat ext , 

2 data( 8) dec(l5,6)\ 

/* LCLINP.TXT - THIS IS PART OF ALL FUNCTIONS AND WAS 

FIXED FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. 

del 

1 lA based ( ia_pt r ) , 

2 I AP float , 

2 lAD float, 

1 IB based( ib_ptr ) , 

2 IBP float, 

2 IBD float, 

1 IC based( ic_ptr ) , 

2 ICP float, 

2 ICE float, 

1 ID based( id ptr), 

2 IDP floatT 
2 lED float , 

1 IE based ( ie_pt r ) , 

2 lEP float , 

2 lED float; 

/* INPUT.AAE */ 
input(l)=UN \ 

input(2) = TWO \ 
input(3) = THREE \ 
input (4) = FOUR \ 
input(5) = FIVE \ 

/* FUNC.AAE -/ 

AD(1 ) = IAP 

/- OUT.AAE */ 

deta( 4) = f tc(Al( 1 ) )\ 

put skip list('data( 4) = ',AD(1))\ 



Figure 6E 

Function Declarations and Equations 
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/* rCLFUN.BAE */ 



AI2 ;proc(code) returns (ptr )\ 

del 

nr_inputs fixed(7) static init(l), 

X2N1 ext entry(f ixed(7) ) returns ( ptr ) , 

TWO ext entry ( fixed (7) ) returns (ptr ) , 

THREE ext ent ry( f ixed ( 7 ) ) returns (ptr ) , 

FOUR ext entry( f ixed (7) ) returns ( ptr ) , 

FIVE ext entry( f ixed(7) ) re turns ( pt r ) , 

cod(5) fixed(7) static in i t (2 , 1 , 1 » 1 , 1 ) , 

F float static init( 0.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data( 8) dec(15,6)\ 

/* LCLINP.TXT - THIS IS PART OF ALL FUNCTIONS ANE WAS 

FIXED FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. -/ 

del 

1 lA based (ia_ptr ) , 

2 lAP float , 

2 lAE float, 

1 IB based ( ib ptr ) , 

2 IBP floatT 
2 I2E float, 

1 IC based ( ic ptr), 

2 IC? float; 

2 ICE float, 

1 IB based (id ptr), 

2 lEP float; 

2 lEE float, 

1 IE based ( ie_ptr ) , 

2 lEP float, 

2 ISE float; 

/- IN PUT. BAD -/ 
input(l) = X2N1 \ 
input ( 2 ) = TWO \ 

input(3) = THREE \ 
input (4) = FOUR \ 
input (5) = FIVE \ 

/* FUN C. BAD */ 

AE12)=IAE 

/* OUT. BA E -/ 

data( 5) = ftc(AD(2))\ 

put Skip list('data( 5) = ',AI(2))\ 



Figure SE 

Func tion Declarations and Equations 



es 



/* ECLFUN.Air */ 

EK :proc(code) returns (ptr)\ 

del 

nr_inputs fiied(7) static init(0), 

ONE ext entry( f iied (7 ) ) returns (ptr ) , 

TWO ext entry( f ixed(7) ) returns (p tr ) , 

THREE ext ent ry ( f ixed ( 7 ) ) returns (ptr ) , 

FOUR ext entry( f ixed (7 ) ) returns (ptr ) , 

FIVE ext entry(f ixed(7) ) returns (ptr ) , 

coa(5) fixed(7) static ini t ( 1 ,1 , 1 , 1 , 1 ) , 

F float static init( 0.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data( 3) dec(15,6)\ 

DCLINP.TXT - THIS IS PART OF ALL FUNCTIONS ANE WAS 
FIXEE FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. */ 

del 

1 lA based ( ia_ntr ) , 

2 lAP float, 

2 lAE float , 

1 13 baseddb ptr), 

2 IBP floatT 
2 IBE float, 

1 IC based( ic ptr ) , 

2 ICP floatT 
2 ICE float, 

1 IE based (id ptr ) , 

2 lEP floatT 
2 ILL float, 

1 IE based ( ie_ptr ) , 

2 lEP float, 

2 lEE float; 

/* INPUT. AEE */ 
input(l) = ONE \ 

input(2) = TWO \ 

input (3) = THREE \ 
input(4) = FOUR \ 
input(5) = FIVE \ 

FUNC.AEE */ 

?=(AD(1)-AE(2) )-CF 

/* OUT. AEE */ 

data( 6) = ftc(F )\ 

put skip list('data( 6) = ' ,1 )\ 



Figure SF 

Function Eeclarations ana Equations 
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/* ECLFUN.BEE */ 

E2K :proc(code) returns (ptr )\ 

del 

nr_inputs fiied(7) static init(l), 

EK ext entry( f ixed(7) ) returns (ptr ) , 

TWO ext ent ry ( f ixed (7 ) ) returns i,ptr ) , 

THREE ext entry( f ixed( 7 ) ) returns ( ptr ) , 

FOUR ext entry(f ixed(7) ) returns (ptr ) , 

FIVE ext entry(f ixed(7) ) returns ( ptr ) , 
cod(5) fixed(7) static init ( 1 ,1 , 1 , 1 , 1 ) , 

F float static init( 0.0000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data( 8) dec(15,6)\ 

ECLIMP.TXT - THIS IS PART OF ALL FUNCTIOMS AND WAS 
FIXEE FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. -/ 

del 

1 lA based ( ia_pt r ) , 

2 lAP float , 

2 lAD float, 

1 IB based ( ib_ptr ) , 

2 IBP float, 

2 IBD float, 

1 IC based(ic_r)tr), 

2 ICP float,' 

2 ICD float , 

1 IE based ( id_ptr ) , 

2 lEP float, 

2 IDE float, 

1 IE based ( ie_pt r ) , 

2 lEP f 1 oat , 

2 lED float; 

/- INPUT. BEE */ 
input (1) = EK \ 

input (2) = TWO \ 
input (3) = THREE \ 
input(4) = FOUR \ 
input(5) = FIVE \ 

/* FUN C. BED */ 

F=(IAP-IAD)/(T-T )-.5-0 
/* OUT. BEE 

data( 7) = ftc(F )\ 

put ski? list('data( 7) = ',F )\ 



Figure 3G 

Function Declarations and Equations 
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3 






/- ECLFUN-AEA -' / 

lAl ;proc(code) returns (ptr )\ 

del 

nr_inputs fixed(7) static init(l), 

E2K ext ent ry( f ixed (7 ) ) returns (ptr ) , 

Ta'O ext entry (fixed (7 ) ) returns (ptr ) , 

TEHEE ext entry(f ixed(7) ) returns (ptr ) , 

FOUR ext entry ( fixed (7) ) returns (ptr ) , 

FIVE ext ent ry (f ixed (7 ) ) returns (ptr ) , 

cod(5) fiied(7) static in i t ( 1 , 1 , 1 , 1 , 1 ) , 

F float static init( 0.2000), 

0 float static init( 0.0000)\ 
del 

1 dat ext, 

2 data( 3) dec(l5,6)\ 

ECLINP.TXT - THIS IS PART OF ALL FUMCTIONS ANL WAS 
PIXEL FOR THE THESIS CASE OF NUMERIC 
CONTROL SYSTEMS. -/ 

id 

1 I A based ( ia_pt r ) , 

2 lAP float, 

2 lAD float , 

1 IB based ( ib_ptr ) , 

2 IB? float, 

2 IBD float , 

1 IC based (ic_ptr) , 

2 I CP float, 

2 ICC float, 

1 II based (id ptr), 

2 IIP floatT 
2 IBD float , 

1 IE based ( ie ptr), 

2 ISP floatT 
2 lED float; 

/- INPUT. AT A -/ 
input(l) = E2K \ 
input (2) = TWO \ 

input (3) = THREE \ 
input (4) = FOUR \ 
input(5) = FI'i/3 \ 

/- FUNC.ALA -/ 

LA(1)=IAP 

/- OUT. ABA -'V 

iata( S) = ftc(BA(l))\ 

put skip list('data( 8) = ',LA(1))\ 



Figure £H 

Function Leclarations and Equations 
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/* LATAINJO.TXT -/ 

timel fixed static init( 100), 
nr_data fixed static init( 3), 
nr_aa fixed static init( 3), 
nr_ad fixed static init( 2), 
nr_dd fixed static init( 2), 
nr~da fixed static init( 1), 
name( 5) char(S) static initl 
'UN ' , 

'XIN ' , 

'X2N1 ' , 

'AEl ' , 

'AB2 ' , 

'EK ' , 

'E2K ' , 

'CAl ' ) 



Figure 9 



EAT A INFO .TXT 



at *lr :• 

* » ^ I • • . 1 ; 

■ U! ■ «• .-si.ff':- 

-^.^0 '• S *! 






I M , I t 



/- FORM. SUB -V 
XSUE 

REFORM AA.TXT AA .TTT 
REFORM AE.TXT AD. TTT 
REFORM EE. TXT EE. TTT 
REFORM EA.TXT EA .TTT 
PLI MONITOR 



REFORM INPUT. AAA INPUT. TTT 
REFORM ECLFUN.AAA ECLFUN.TTT 
PIP FUNG. TXT=FUNC .AAA 
REFORM OUT. AAA OUT. TTT 
PLI AA 

PIP UN .REL=AA.REL 

REFORM INPUT. BAA INPUT. TTT 
REFORM ECLFUN.BAA ECLFUN.TTT 
PIP FUNG. TXT=FUNC .BAA 
REFORM OUT. BAA OUT. TTT 
PLI AA 

PI? X1N.RSL=AA.REL 

REFORM INPUT. CAA INPUT. TTT 
REFORM ECLFUN.CAA ECLFUN.TTT 
PIP FUNG. TXT=FUNC .CAA 
REFORM OUT. CAA CUT. TTT 
PLI AA 

PIP X2N1. REL=AA .REL 

REFORM INPUT. AAD INPUT. TTT 
REFORM ECLFUN.AAE ECLFUN.TTT 
PIP FUNC.TXT=FUNC.AAE 
REFORM OUT.AAE OUT. TTT 
PLI AE 

PIP AE1.REL=AE.REL 

REFORM INPUT. BAE INPUT. TTT 
REFORM ECLFUN.BAE ECLFUN.TTT 
PIP FUNC.TXT = FUNC .BAD 
REFORM OUT. BAE OUT. TTT 
PLI AD 

PI? AE2 .REL=AL. HEL 



Figure 12A 



FORM .3'JB 



74 



REFORM input. ATE INPUT. TTT 
REFORM DCLFUN.ADD DCLFUN.TTT 
PIP FUNG .TXT=FUNC .AEE 
REFORM OUT .AEE OUT. TTT 
PLI DD 

PIP EK.REL = EE .REL 

REFORM INPUT. BEE INPUT. TTT 
REFORM ECLFUN.BEE ECLFUN.TTT 
PIP FUNG. TXT=FUNC. BEE 
REFORM OUT. BED CUT. TTT 
PLI EE 

PIP E2E.REL=EE.REL 

REFORM INPUT .AEA INPUT. TTT 
REFORM EGLFUN.AEA ECLFUN.TTT 
PIP FUNG. TXT=FUNG. AEA 
REFORM OUT. AEA OUT. TTT 
PLI EA 

PIP DAI .RSL=DA .REL 



PLI PRINT 
LINK PRINT 

LINK ANSjv'SR=MONITOR ,UN ,X1N ,X2N1 ,AD1 ,AD2 ,EA1 , 



Figure 103 
FORM. SUB 
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E. TEST RESULTS ANE PROGRAM EXECUTION 

As mentioned previously, once the xzwgzau ANS^VER.COM is 
developed, all that is necessary is to enter ANS-i/SF. to the 
CP/M operating system promnt. Figure 4B and 9 should he 
referenced during this discussion. Figure 11 is a partial 
printing of what will he on the display when the urogram is 
executing. The printout is for the sample problem but with 
the EELAY set to zero and the ELIMIT shortened. The first 
thing to note is that the functions are not executed ir f rder 
and some are executed several times during each clock. Time 0 
actually runs from ECLOCX = O to ACLCCK = 1. By associating 
the data element number with the order of the function names 
in Figure 9 the user can determine which function is oeing 
calculated or called. As an example the first element 
returning an answer in Figure llA is element 4 which equates 
to All in Figure 9. The next elements are 5, 6, and 7 that 
equate to AD2, EX, and S2X respectively. An important note is 
that since there is no delay, function 3 (EAl) or the di;Cital 
to analog function is calculateu or called at time 0. This 
will not be the case when discussed later on. The only 
function that has a value is element 1 or JN. There is no 
further changes in the data until DCLCCS = 1 which is the 
same time as ACLOCK = 0 for the second time. At this point 
the analog to digital conversion has taken place for the 
first time, the computer does internal calculation, dieital 
to analog is complete since there is no delay, and the analog 
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plant starts to integrate the signal. This is really the 
start of the sample since the computer is operating on the 
initial conditions for the first ten counts of the analog 
clock. By the first sample period after the initial values 
(rCLOCK = Z) the controller output changes to -1 and the 
plant output is at the half way value of the input. By the 
second sample period after the initial conditions 'ECLOCK = 
3) the plant output is equal to the forcing function i-vnut 
and the controller output is 0. It should he noted that in 
all programs the zero point of the sample period always 
starts at ECLOCS = 0 so that the initial conditions ca-’ De 
graphed. The above example is a theoretical non realiralle 
simulation of a control system with no calculation delays. 
The results obtained are identical to classical texiiock 
exampl es . 

Figure 12 is the identical problem discussed with the 
delay equal to zero but it shows what the display file looks 
like after the utility program PHINT.FLI is run. The results 
in Figure 12 are obviously more readable sinc<= the function 
name and final value are not presented until after all 
calculations are complete. 

Figure 13 is a partial printing of the results wr.en using 
the full example presented above with a forty percent delay. 
This means there is four ACLOCK counts after the ECLCCi 
before the digital to analog values are available. The 
important points to note are that the digital to analog 
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values (DAI or function 8) do no change until the analog 
clock ( ACLOCK ) is equal to 4 during each digital clock, the 
overshoot is only about 2.5 percent, and it takes five times 
as long to reach stability or the input signal. These results 
are much more representative of a real world problem ana the 
problems that would be encountered with a true hardware 
impl ementat ion . 
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ANSWER. C0(^. Runtime Results 
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Figure 12A 

PRINT. PLI With Zero Eel ay 
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PRINT. ?LI Xith Zero Delay 
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PRINT. PLI With Forty Percent Eelay 
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Figure 12B 



PRINT. PLI 'A'iih Fortj'’ Percent Eela/ 
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Figure 13C 
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PLI th Forty Percent Delay 
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X2N1 


0 . 999450 


ADI 


1.000000 


AE2 


1.000391 


EK 


-0 .000262 


E2E 


0 .001065 


EAl 


0.001068 


TIME=100 


UN 


1 .000000 


XIN 


-0 .000139 


X2N1 


0.999431 


ADI 


1.000000 


AE2 


0.999450 


EX 


0 .000368 


S2K 


0 .000096 


EAl 


3.001068 


PRINT. PLI 


Figure 13E 
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VIII. CONCLUSIOMS 



A major point of this thesis was to demonstrate the 
feasibility of writine a computer program that can be used on 
a microcomputer by a novice programmer to develop cor.trol 
system models, simulation, and realization. This type ci 
system was not possible in the past due to system memory 
limitations and. unacceptable software programs. The approach 
taken in this project was to use commercially availaole hi£h 
level programming languages to reduce the coding time. As vas 
discussed in this thesis, the languages are not sufficient in 
themselves and more than one lanmace was needed to ccmnlete 
the project. It is apparent that a language that has the 
features of a good relational database couple-l with extremely 
good mathematical functions and the ability to communi'ate 
with the computer hardware is a necessity for this type r: 
system to be a useful tool. 

The importance of recognizing that control systems can be 
modeled and manipulated in a database cannot be ov-riooice:. 
ly carrying this approacn to the limit it is possible to 
develop a large database of control functions that could be 
recalled and used like a library. Once tue database is formed 
it could be manipulated by artificial intelligence. T'::e 
computer itself could be told to develop a control system to 
match a desired set of waveforms using the information stereo 



or at least find a best fit solution, f^issing parameters 
could be searched for, or best fit solutions ccula possibly 
be provided as easily as using the the present database 
systems to find phone numbers and names. 

It is significant that source code was used for the 
intermediate language. A special purpose program coula have 
been written that provided direct function calls to specific 
nardware; however, by using high level language sour'e ' ode 
the program will run on many different machines just by usinf- 
their * resident compilers. This was demonstrated by 
transporting the prog’^ams to four different me 'hir.es with 
different operating environments. Ey using source code the 
problem was kept at a high level of abstraction. The orogramis 
the computer needed* to write were similar to English which 
helped lessen the confusion while developing more than o.na 
hundred files and programs. 

fact that became very apparent early in the urogram 
development was the need for a correction factor to ta-'e into 
account the delays in the system due to calculations and 
conversions. This problem is aadressed in many models but it 
is seldom considered in most computer solutions sin'^e the 
computation time is usually insignificant. It was neterminei 
that solving for all possible delays was next to impossible 
especially when implemented since the hardware itself would 
be variable. A correction factor is provided that works very 
well with the example problem. Eurt.ner investigation needs to 



be done using other functions to test the results; however, 
it is felt that the factor will work with any control system 
designed to have a zero steady state error signal. The basic 
approach is to have the designer develop the model assuming 
there are no delays and then estimate the delays to be 
applied to the correction factor. If the factor does not work 
then try to develop one from other experience. If this fails 
then use benchside engineering, make estimates and watch what 
the results do. 

Appendix C contains several graphs of data covering the 
example problem. A variety of delays were introduced an: the 
results with and without the correction factor are proviiied. 
The problefn was initially designed assuming no delay and tnen 
corrected for with the correction factor. It should be noted 
that in all cases, even with delays up to £9 pe-rent, 
acceptable results were produced within five times the 
initial sample period design. This is significant since the 
uncorrected functions had severe problems even with the ‘■mall 
delays. When using the correction factor the longer delays 
corrected better and it is felt that this is due ^'o 
outputting a control signal closer to the time the input was 
going to sample the result of the control. If a signal is 
output and sampled immediately, the probability that tie 
response to that signal is the one observed is highei than if 
a long delay were encountered possibly resulting in severe 
data error. Sven though the system was designed for a step 

S£ 



input there is evidence that the correction factor works very 
well for a ramp input. This is also shown in Appendix C. 

As a final note the power of a granhical dynamic display 
for the designer is invaluable. Reading through many pages of 
data makes it very difficult to see where the significant 
changes are being made. A real time plot will cut design time 
by an order of magnitude. 



APPENEIX A 



E3ASS II PROGRAMS 



CINIT EBF 
CMENU CME 
CONTROL CME 
DBASEAPP OVH 
EBASSMOE OVH 
EBASESRT OVR 
PARINEEX NEX 
FORMl CME 
MENU FMT 
STRT MEM 
XSUE COM 
PROC CME 



CINIT FMT 
CMENU FMT 
PARAMS EBF 
EBASEBRO OVR 
EBASSMSC OVR 
EBASSTTL OVR 
SFORM SUB 
INIT EBF 
MENU ZIP 
SUBMITl COM 
SHORT CME 
EATA CME 



CINIT ZIP 
CMENU ZIP 
EATAl CME 
EBASSJOI OVR 
EEA3EM3G COM 
LBASEUPE OVR 
ECLFUN CME 
INIT MEM 
MENU ZPR 
ZIP COM 
ECLM.ON CME 
FORM CME 



CINIT ZPR 
CMENU ZPR 
EBASE COM 
EBASEMAI CVR 
EEASERPG OVR 
EBINSTL COM 
FORM SUB 
MENU CME 
STRT EIF 
ZIPIN COM 
MONECL CME 
RSFCRM CME 



Figure 14 

Programs Required for Eatabase Operation 
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=!* MENU. CMC 
SET TALK OFF 

RESTORE FROM INIT See Figure 2 

USE INIT See Figure 3 

SET FORMAT TO MENU See Figure 15A and 15E 

EO WHILE T 
READ 
DO CASE 

CASE ACTION = ' 0 ' 

ACCEPT 'DO YOU REALLY WANT TO QUIT? Y/N ' TO CN 
IF CN = 'y' 

CUIT 

SNDIF 

CASE ACTION = 'l' 

RESTORE FROM INIT 
USE INIT 
COPY TO PARAMS 
USE PARAMS 
DO CMENU 

CASE ACTION = '2' 

ACCEPT 'ENTER NAME OF OLD PROGRAM FILE ' TO PARAM 

RESTORE FROM 'LPARAM 

USE S.PARAM 

COPY TO PARAMS 

USE PARAMS 

EO CMENU 

CASE ACTION = '3' 

ACCEPT 'ENTER NAME OF FILE TO SAVE PARAMAT2R3 IN ' TO PARAM 
SAVE TO 'iPARAM 
COPY TO &PARAM 
CASE ACTION = '4' 

EISP ALL FUNCNAME 

ACCEPT 'PRESS RETURN TO CONTINUE ' TO CN 
CASE ACTION = '5' 

ACCEPT 'ENTER NAME OF OLE PROGRAM FILS ' TO PARAM 
RESTORE FROM SPARAM 
USE SPARAM 
COPY TO PARAMS 
USE PARAMS 
DO SHORT 
ENDCASE 
ENEEO 



c 



file MENU 



JU «v 

^ ^ 't* 



ENTER 0 TO QUIT 

ENTER 1 TO START A NEW PROBLEM 
ENTER 2 TO RECALL AN OLE PROBLEM 
ENTER 3 TO SAVE PRESENT PARAMATExRS 
ENTER 4 TO DISPLAY PRESENT EUNCTIONS 
ENTER 5 TO CHANGE INITIAL CONDITIONS ONLY 
( This will allow the pro.gratn 
to compile much faster ) 

ENTER ACTION #ACTION 



Figure ISA 



MENU 





MENU 


FMT 








0 


7,17 


SAY 


"enter 


0 TO 


quit" 


0 


8,17 


SAY 


"enter 


1 TO 


START A NSW PROBLEM" 


0 


9,17 


SAY 


"enter 


2 TO 


RECALL AN OLD PROBLEM" 


0 


10,17 


SAY 


"enter 


3 TO 


SAVE PRESENT PARAMATERS" 


0 


11,17 


SAY 


"enter 


4 TO 


DISPLAY PRESENT FUNCTIONS 


0 


12,17 


SAY 


"enter 


5 TO 


CHANGE INITIAL CONDITIONS 


0 


13,28 


SAY 


"( This 


will 


allow the program" 


0 


14,30 


SAY 


to compile 


much faster )" 


0 


16,23 


SAY 


ENTER 


ACTION 


0 


16,36 


GET 


ACTION 







Figure 15B 
MENU Format 



* CMEMU.CME 

SET FORr^AT TO CMSNU * See figure 16A and 163 
EO WHILE T 
RSAE 
EC CASE 

CASE ACTION = '0' 

PACK 

INDEX ON FUNCTYPE TO PAfilNEEX 
USE PARAMS INEEX PARINEEX 
DO CONTROL 
CASE ACTION = 'l' 



ERASE 

LISP ALL FUNCNAME 
ACCEPT 'ENTER RECORD 


NUMBER 


TO SLIT 


' TO CN 


EEIT &CN 

CASS ACTION = '2' 
APPEND FROM INIT 


- See fi 


gure 3 




EEIT # 

CASS ACTION = '3' 

ERASE 

EI3P ALL FUNCNAME 
ACCEPT 'ENTER RECORD 


NUMBER 


TO DELETE 


' TO CN 


DELETE RECORD ECN 
CASE ACTION = '4' 

ERASE 

EI3P ALL FUNCNAME 
ACCEPT 'ENTER RECORD 


NUMBER 


TO RECALL 


' TO CN 


RECALL RECORD S.GN 
CASE ACTION = 'S' 

SET FORMAT TO CINIT 


* S ee 


figure 1?A 


and 17B 


RSAE 

SET FORMAT TO CMENU 


- See 


figure 16A 


and 16E 


ACCEPT 'PRESS RETURN 


TO CONTINUE ' TO 


CN 


CASS ACTION = '6' 

SET FORMAT TO MENU 


* See 


figure 15A 


and 151 


RETURN 







ENECASS 

ENEDO 

RETURN 
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file CMENU 



ENTER 

ENTER 

ENTER 

ENTER 

ENTER 

ENTER 

ENTER 



- CMENU.FMT 
0 3,15 3aY 
0 4,15 SAY 
0 5,15 SAY 
0 6,15 SAY 
0 7,15 SAY 
0 8,15 SAY 
0 9,15 Say 
0 12,22 SAY 
0 12,35 GST 



0 TO COMPILE A PROGRAM 

1 TO ELIT TRANSFER FUNCTIONS 

2 TO AIL A TRANSFER FUNCTION 

3 TO DELETE A TRANSFER FUNCTION 

4 TO RECALL A DELETED TRANSFER FUNCTION 

5 TO CHANGE INITIAL CONDITIONS 

6 TO RETURN TO THE MAIN PROGRAM 



ENTER ACTION #ACTION 



Figure 16A 
CMSNU 



enter 

■’enter 

"enter 

"enter 

"enter 

"enter 

"enter 

"enter 

ACTION 



0 TO COMPILE A PROGRAM 

1 TO EDIT TRANSFER FUNCTIONS" 

2 TO ADD A TRANSFER FUNCTION" 

3 TO DELETE A TRANSFER FUNCTION" 

4 TO RECALL A DELETED TRANSFER FUNCTION" 

5 TO CHANGE INITIAL CONDITIONS" 

6 TO RETURN TO THE MAIN PROGRAM" 

action" 



Figure 16B 
CMENU Format 
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Jile CINIT ** 



DLIMIT #DLIMIT (This is the number of times the computer 

will do calculations, ie. The number of 
sample periods . ) 

ALIMIT #ALIMIT (This is the number of times each analog 

calculation will be done during each sample 
period. Effects how accurate the analog plant 
is modeled . ) 

IELAY #EELAY (This is the number of analog calculations done 

before the digital to analog is availible. ) 

T #T (This is the sample rate of the digital system. 

You must provide it if you are going to use it 
in your equations.) 

T1 #T1 (This is the sample rate of the analog system. 

It is usually the same as T above during design; 
however, for simulation it usually works out to 
be T divided by ALIMIT since you usually want a 
more accurate plot of the analog nlant to see 
what happens between the digital sample periods. 
You must provide T1 if you use it in your equa- 
tions.) (This is effectively double rate sampling. 



Figure 17A 
CINIT 



G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 



CINIT, 


,FMT 


2 


. 3 


SAY 


2 


,10 


(JET 


2 


,24 


SAY 


3 


,25 


SAY 


4 


,25 


SAY 


5 


, 3 


SAY 


5 


,10 


GET 


5 


,24 


SAY 


6 


,25 


SAY 


7 


,25 


SAY 


8 


,25 


SAY 


9 


, 3 


SAY 


9 


,10 


GET 


9 


,24 


SAY 


10 


,25 


SAY 


11 


, 4 


SAY 


11 


,10 


GET 


11 


,24 


SAY 


12 


,25 


SAY 


13 


,25 


SAY 


14 


, 4 


SAY 


14 


,10 


GET 


14 


,24 


SAY 


15 


,25 


SAY 


16 


,25 


SAY 


17 


,25 


SAY 


18 


,25 


SAY 


19 


,25 


SAY 


20 


,25 


SAY 


21 


,25 


SAY 



"dlimit" 

ELIMIT 

"(This is the number of times the computer" 

"will do calculations, ie. The number of" 

"sample periods . )" 

"alimit" 

ALIMIT 

’](This is the number of times each analog" 
"calculation will be done during each sample" 
"period. Effects how accurate the analog plant" 

"is modeled . ) " 

"eelat" 

DELAY 

]|(This is the number of analog calculations done" 
"before the digital to analog is availible . 1 ’’ 

T 

T 

"(This is the sample rate of the digital system." 
"You must provide it,,if you are going to use it" 

"in your equations.)" 

"ti" 

T1 

"(This is the sample rate of the analog system." 

*[lt is usually the same as T above during design;" 
"however, for simulation it usually works out to]’ 
]]be T divided by ALIMIT since you usually want_ a" 
]]more accurate plot of the analog plant to see" 
"what happens between the digital sample periods." 
]]You must provide Tl if you use it in your equa-" 
"tions . ) (This is effectively double rate sampling.) 



Figure 17B 
CINIT Format 
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* CONTROL. CM! 

EO ECLMON 

STORE 0 TO INI 

LOCATE FOE FUNCTYPE = 'AA' 

EO PROC 

STORE 0 TO INE 

LOCATE FOR FUNCTYPE = 'AE' 

EO PROC 

STORE 0 TO INE 

LOCATE FOR FUNCTYPE = 'IE' 

EO PROC 

STORE 0 TO INE 

LOCATE FOR FUNCTYPE = 'EA' 

EO PROC 

EO EATA 

EO FORM 

QUIT TO 'SUBMITI FORM' 

SET TALK ON 
RETURN 
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* DCLMON.CMD 

SET ALTERNATE TO B:ECLt^ON - B: not needed if single disk 

SET ALTERNATE ON 

LOCATE 

EO WHILE .NOT. EOF 

? ' '+FUNCNAME+' ext entry (fixed (7 ) ) returns (ptr ), ' 
CONTINUE 
ENEEO 

COUNT FOR SAVEDATA TO DATAl 

7 ' T float static ext ini t { '+s t r ( T ,8 ,4) + ' ) , ' 

? ' T1 float static ext init( ' + str (T1 ,8,4)+ ' ) , ' 

? ' ELIMIT fixed static init('+ str (BLIMIT ,4) + ' ) , ' 

? ' ALIMIT float static init('+ str(ALIMIT ,4) + ' ) , ' 

? ' delay float static init('+ str(EELAY , 4)+ ' 

7 ' 1 dat ext , ' 

? ' 2 data( '+str (DATAl ,3) + ') dec(15,6)' 

SET ALTERNATE OFF 
STORE 'AA' to T:NAME 
DO MONDCL 

STORE 'AD' TO TrNAME 
EO MONDCL 

STORE 'EE' TO T;NAME 
EO MONDCL 

STORE 'EA' TO T:NAME 

EO MONECL 

RETURN 



* MONECL.CMD 

SET ALTERNATE TO B:&T:NAMS 3: not needed if single disk 
SET ALTERNATE ON 
LOCATE FOR FUNCTYPE = T:NAME 
IF EOF THEN 
? 

ELSE 

EO WHILE .NOT. EOF 

? 'dummy_ptr = ' + funcname + ' ^code)\ ' 

CONTINUE 

ENDDO 

SNEIF 

SET ALTERNATE OFF 
RETURN 
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* PROC.CME 
EO WHILE .NOT. EOF 
STORE IND+1 TO INE 

STORE 'ECLFUN. '+$(FUN,INE ,1)+FUNCTYPE TO T:NAME 

SET ALTERNATE TO B:S<T:NAME * P: not needed if single disk 

SET ALTERNATE ON 

? FUNCNAME+ ' :proc ( code) r eturns ( ptr ) \ ' 

7 'del' 

? 'nr_inputs fixed(7) static init ( '+HR : INPUTS + ' ) , ' 

? INPUT:A+' eit entry (f ixed (7 ) ) returns (ptr ), ' 

? INPUT:B+' ext en try ( f ixed (7 ) ) returns ( ptr ), ' 

? INPUT:C+' ext entry (fixed (7 ) ) returns (ptr ), ' 

? INPUT;E+' ext entry(f ixed(7) ) returns(ptr) , ' 

? INPUT:E+' ext entry(f ixed (7) ) returns (ptr ), ' 

? 'cod(5) fixed(7) static init ( ' + /UP2E + ' , '+B1P2E + ; 

', ' + C1P2E+', ' + ElP2E+',' + ElP2E+') , ' 

STORE STR(F5,8,4) TO TEMPOUT 
? '? float static ini t ( '+TSKPOUT+ ' ) , ' 

STORE STR(FKE,8.4) TO TEMPOUT 
7 '0 float static ini t ( '+TEMPOUT+ ' )\ ' 

7 'del' 

7 ' 1 dat ext , ' 

7 ' 2 data( '+STfi(datal ,3) + ') dec(15,6)\' 

SET ALTERNATE OFF 

STORE 'INPUT. '+$(FUN, INE, 1)+FUNCTYPE TO T:NAME 
SET ALTERNATE TO B:SiT:NAME B: not needed if single disk 
SET ALTERNATE ON 
7 'input(l) = '+INPUT:A + '\ ' 

7 'input(2) = '+INPUT :B+'\' 

? 'input(3) = '+INPUT:C+'\' 

? 'input(4) = '+INPUT:E+'\' 

7 'input(5) = '+INPUT:£+'\' 

SET ALTERNATE OFF 

STORE 'FUNG .'+$ (FUN , INE ,1 )+FUNCTYPE TO TrNAME 

SET ALTERNATE TO BtSTtNAMS B: not needed if single disk 

SET ALTERNATE ON 

7 EQUATION 

SET ALTERNATE OFF 

STORE 'OUT.'+$(?UN,rNE,l)+FUNCTYPE TO TrNAME 
SET ALTERNATE TO BrETrNAMS B: not needed if single disk 
SET ALTERNATE ON 
IF SAVEEATA 

STORE ELEMSNTS+1 TO ELEMENTS 

? 'data( ' + STR(ELEMENTS ,3) + ') = f tc ( '+EATANAMS+ ')\ ' 

7 [put skip list( 'aata(]+STR(ELEMENTS ,3)+ [) = ',] 

?? EATANAME+ ')\' 

ELSE 

7 

ENEIF 

SET ALTERNATE OFF 
CONTINUE 
ENEEO 
RETURN 
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DATA.CMD 

SET ALTERNATE TO EtEATAINFO * B: not needed if single disk 
SET ALTERNATE ON 

7 'timei fixed static init( ' + STR(ELIMIT*ALIMIT,4) + ') , ' 

? 'nr data fixed static ini t ( '+STR( EATAl ,3 ) + ' ) , ' 

C0UNT"FCR SAVEDATA .and. (FUNCTYPE = 'kk') TO INL 
7 'nr_aa fixed static ini t ( '+STR ( I ND, 2 )+ ' ) , ' 

COUNT FOR SAVEDATA .AND. (FUNCTYPE = 'AD') TO IND 
7 'nr ad fixed static ini t ( '+STR( IND , 2 )+ ' ) , ' 

C0UNT"F0R SAVEDATA .AND. (FUNCTYPE = 'DD') TO IND 
7 'nr dd fixed static ini t ( '+STR( I ND ,2 )+ ') , ' 

COUNT'FOR SAVEDATA .AND. (FUNCTYPE = 'DA') TO IND 
7 'nr_da fixed static ini t ( '+STR ( IND ,2 )+ ' ) , ' 

7 'name( '+STR(DATA1 ,3) + ') char(£) static init(' 

STORE 'AA' TO T : NAME 
DO DATAl 

LOCATE FOE FUNCTYPE = 'AA ' 

IF .NOT. EOF 
77 ' 

ENDIF 

STORE 'AD' TO T:NAME 
DO DATAl 

LOCATE FOR FUNCTYPE = 'AD ' 

IF .NOT. EOF 

?? '/ 

ENDIF 

STORE 'DD' TO TtNAME 
DO DATAl 

LOCATE FOE FUNCTYPE = 'DD ' 

IF .NOT. EOF 
77 ',' 

ENDIF 

STORE 'DA' TO TtNAME 
DO DATAl 
77 ')' 

SET ALTERNATE OFF 
RETURN 



* DATA1.CMD 

LOCATE FOR SAVEDATA .AND. (FUNCTYPE = TtNAME) 
DO WaiLS .NOT. EOF 
? C']+FUNCNAME + ['] 

CONTINUE 
IF .NOT. EOF 
7? ',' 

ENDIF 

ENDDO 

RETURN 
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* FORM.CME 

SET ALTERNATE TO FORM. SUB 
SET ALTERNATE ON 

? ' 3 :' * B: not needed if single disk 

? 'XSUB' 

? 'REFORM AA.TXT AA.TTT' 

? 'REFORM AB.TXT AB.TTT' 

? 'REFORM DD.TXT DD.TTT' 

? 'REFORM EA.TXT EA.TTT' 

? 'PLI MONITOR' 

STORE 'AA' to T;NAME 
EO ECLFUN 

STORE 'AE' to T;NAME 
DO DCLFUN 

STORE 'EE' TO T:NAME 
EO ECLFUN 

STORE 'DA ' TO T:NAME 
EO ECLFUN 
? 'PLI PRINT' 

? 'LINK PRINT' 

? 'LINK ANSWER=MONlTOR' 

LOCATE 

DO WHILE .NOT. EOF 
?? ', '+TRIM(FUNCNAME) 

CONTINUE 

ENDDO 

? 

SET ALTERNATE OF? 

RETURN 



* DCLFUN.CMD 

LOCATE FOR FUNCTYPE = T:NAME 
STORE 0 TO INE 
EO WHILE .NOT. EOF 
STORE IND^l TO IND 

STORE $ (FUN ,INE, 1 )+?UNCTYPE TO EXT:NAME 
EO REFORM 
CONTINUE 
ENEEO 
RETURN 



* RSFORM.CMD 

? 'REFORM INPUT . '+SXT :NAME+ ' INPUT.TTT' 

? 'REFORM ECLFUN. '+EXT :NAME+' ECLFUN .TTT ' 

? 'PIP FUNG .TXT=?UNC . '+EXT: NAME 
? 'REFORM OUT. '+EXT :NAME+' OUT .TTT ' 

? 'PLI '+FUNCTYPE 

? 'PIP '+TRIM(FUNCNAME)+'.RSL='+?UNCTYFS+'.REL' 
RETURN 
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* SHORT.CMD 



SET FORI^AT TO CINIT * See figure 17A and 17B 
READ 

SET TALK OFF 

SET ALTERNATE TO BcLCLMON *B : not needed if using a sinsle disk 

SET ALTERNATE ON 

LOCATE 

10 WHILE .NOT. EOF 

? ' '’+FUNCNAME+ ' ext en try (f ixed (7 ) ) returns ( ptr ), ' 

CONTINUE 

ENLLO 

COUNT FOR SAVEEATA TO EATAl 

? ' T float Static ext init( '+str (T,8,4) + ') . ' 

? ' T1 float static ext init ( '+str (Tl,8,4)+ ' 

? ' ELIMIT fixed static init('+ s tr( LLIMIT , 4)+ ' ) , ' 

? ' ALIMIT float static init('+ str (ALIMIT ,4) + ' ) , ' 

? ' delay float static init('+ str (BELAY, 4)+ ' 

? ' 1 dat ext,' 

? ' 2 data( '+str(BATAl ,3)+' ) dec(l5,6)' 

SET ALTERNATE OFF 
SET ALTERNATE TO SFORM.SUB 
SET ALTERNATE ON 
? 'XSUB' 

? 'B:' * E: not needed if using a single disk 

? 'PLI MONITOR' 

? 'PLI PRINT' 

? 'LINK PRINT' 

7 'LINK ANSWER=M0NIT0R' 

LOCATE 

LO WHILE .NOT. EOF 
?? ', ' + TRIM(FUNCNAME) 

CONTINUE 

ENBEO 

? 

SET ALTERNATE OFF 
QUIT TO 'SUBMITl 3F0RM' 

RETURN 
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APPENEIX B 



PL/1-S0 PROGRAMS 



A A 


PLI 


DD 


PLI 


PLI0 


OVL 


PIP 


GOM 


REFORM 


SYM 


EE 


TXT 


FUiMC 


AAA 


FUNG 


BAA 


FUNG 


GAA 


FUNG 


AAE 


FUNG 


BAD 


FUNG 


AEE 


FUNG 


BDD 


FUNG 


AEA 


AE 


TTT 


INPUT 


TTT 


XIN 


REL 


AE2 


REL 


EK 


REL 


PRINT 


REL 


ANSWER 


GOM 



AE 


PLI 


LINK 


GOM 


PLIl 


OVL 


REFORM 


PLI 


EGLMON 


TXT 


EA 


TXT 


OUT 


AAA 


OUT 


BAA 


OUT 


GAA 


OUT 


AAE 


OUT 


BAD 


OUT 


AEE 


OUT 


BDD 


OUT 


AEA 


EE 


TTT 


DGLFUN 


TTT 


UN 


REL 


ADI 


REL 


EA 


REL 


PRINT 


GOM 


ANSWER 


SYM 



EA 


PLI 


MONITOR 


PLI 


PLI2 


OVL 


REFORM 


REL 


A A 


TXT 


EGLFUN 


AAA 


DGLFUN 


BAA 


EGLFUN 


GAA 


EGLFUN 


AAE 


EGLFUN 


BAE 


EGLFUN 


ADD 


EGLFUN 


BEE 


EGLFUN 


ADA 


EATAINFO 


TXT 


EA 


TTT 


CUT 


TTT 


X2N1 


REL 


LD 


ReL 


FUNG 


TXT 


PRINT 


SYM 


PRINT 


PLI 



EGLINP 


m y fTi 


PLI 


GCM 


PL I LIB 


IRL 


REFORM 


COM 


A I 


TXT 


I N PUT 


AAA 


INPUT 


BAA 


INPUT 


GAA 


INPUT 


AAD 


INPUT 


BAE 


INPUT 


A. EE 


INPUT 


3IE 


INPUT 


AEA 


A A 


TTT 


MON IT OR 


nEL 


AA 


REL 


AE 


REL 


E2X 


REL 


EAl 


REL 


DATA 


FLT 



Figure 13 

Programs Required for FL/1-60 Operation 
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MONITOR: proc opt ions (main , s tack ( 512 )) ; 
del 

ftc entry (float binary) returns (char ( 17 ) var)» 

del 

TIME fixed ext static init(0), 

DCLOCK fixed ext static init(0), 

ACLOCK fixed ext static init(0), 

AE(8) float ext static ini t ( 0 ,0 ,0 ,0 ,0 ,0 , 2 ,0 ) , 

DA(8) float ext static ini t (0 , 0 , 0 , 0 , 0 ,0 , 0 ,0 ) , 
dummy_ptr ptr, 
code f ixed (7 ) , 

CF float ext, /* correction factor for error signal 
datafile file; 

del 

%include 'dclmon.txt'; 

> 

open f ile(dataf ile ) record sequential output t itle ( 'data .f 1 1 ' ) 
env(b(512) ) ; 

CF = ftc(exp(char(-delay/ALIMIT) ) ); 
do while (DCLOCK <= DLIMIT); 

PUT SKIP LIST(ECLOCK, 'DCLOCK'); /=!' SLIT ^' / 
code = i; 

ACLOCK = 0; 

%include 'ad.ttt'; 
iinclude 'dd.ttt'; 
code = i; 

do while (ACLOCK <= ALIMIT - 1); 

PUT SKIP LIST(ACL0CK, 'ACLOCK', TIME, 'TIME'); /- EEIT - 
if ACLOCK = delay then do; 

^include 'da . ttt ' ; 
end; /* if */ 

%include 'aa .ttt ' ; 

write f ile(dataf ile ) from(dat); 

TIME = TIME + i; 

ACLOCK = ACLOCK + i; 
end; /* do */ 

DCLOCK = DCLOCK + i; 
end; /* do / 
close f ile (dataf i le ) ; 
end monitor; 
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/* AA.PLI 

^include 'dclf un. tt t ' ; 
del 

ftc entry (float binary) returns { cha r( 17 ) var); 
del 

TIME filed ext , 

T float ext, 

T1 float ext, 

input(5) variadle entry (fixed (7) ) re turns ( ptr ) , 
EA fixed statie init(0), 

EAl fixed statie init(0), 

DA(8) float ext, 

ACL0C5 fixed ext, 
eode fixed (7 ) , 
result ptr, 
m(5) label, 

( ia_ptr , ib_ptr , ie ptr, id_ptr, ie_ptr ) ptrj 
iinelude 'delinp . txt"^; 



%include 'input. ttt'; 
if eode = 1 then do; 

if EA < TIME then do; 

if (nr_inputs> 0) S. (nr_inputs< 6) 
then~do ; 



m ( 1 )=la ; 
m(2)=lb; 
m (3 ) = le ; 
m(4:)=ld; 
m ( 5 )=le ; 

goto m ( nr_ inputs ) ; 
end; /* then 
else goto If; 
le: ie_ptr = input (5) 

Id: id_ptr = input (4) 

le: ie_ptr = input (3) 

lb: ib_ptr = input (2) 

la: ia_ptr = input (l) 

If: 0 = f; 

EAl = time; 

%inelude 'fune.txt'; 



(cod(5) ) 
(eod (4) ) 
( eod (3) ) 
( c 0 d ' 2 ) ) 
( cod ( 1 ) ) 



EA = time; 
end; /- if EA 
end; if code then */ 
else 

if EAl < TIME then do; 

0 = f; 

EAl = time; 
end; EAl */ 
^include 'out. ttt'; 
result = addr(F); 
return (result ) ; 
end ; 
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/* AE.PLI ’*'/ 

^include 'dclf un. ttt 



del 

ftc entry (float binary) returns (char{ 17) var); 
del 

EA fixed statie init(0), 

T float ext, 

T1 float ext, 

input(5) variable entry ( fixed (7 ) ) returns ( ptr ) , 
CE f loa t ext , 

AE(8) float ext, 
eode fixed(7), 
result ptr, 
m(5) label, 

(ia_ptr,ib_ptr,ic_ptr,id_ptr,ie_ptr) ptrJ 

^inelude 'dclinp.txt'; 

%inelude 'input. ttt'; 

if EA > 0 then do; 

if (nr_inputs> 0) & (nr_inputs< 6) 
then do; 
tn ( 1 )=la ; 

[Ti(2)=lb; 
m( 3)=le ; 
m(4)=ld; 
m(5)=le; 

goto [Ti( nr_ input s ) ; 
end; then */ 
else goto If; 



le : 


ie_ptr 


= input 


(5) 


( cod ( 5 ) ) ; 


Id: 


id_ptr 


= input 


(4) 


(cod (4 ) ) ; 


le : 


ie_pt r 


= input 


(3) 


(cod (3) ) ; 


lb; 


ib_ptr 


= input 


(2) 


(cod (2) ) ; 


la : 


ia_ptr 


= input 


(1) 


(cod { 1 ) ) ; 


If: 










iinelude 


'f une . txt ' ; 


/- AICO 



end; then */ 
result = addr(F); 

EA = i; 

%inelude 'out.ttt'; 
J 

return (result ) ; 
end ; 



, « e-e-e-s-fj s. i ^_, 

i>ij-ij-ij-ij-ij-ij-ij-ij-ij-ij-^^^^^ ^ y y 1 

3 J;-iJ-iJ-iJ-iJ-iJ-iJ-iJ-iJ-iJ-iJ-iJ-iJ-^^^^^ ]J FI :J rj 



/* DC.PLI */ 

%include 'dclf un . tt t ' ; 
del 

ftc entry (float binary) returns ( char( 17 ) var) 

del 

TIME fixed ext, 

T float ext, 

T1 float ext, 

input(5) variable entry( f ixed (7 ) ) returns(ptr) 
EE fixed static init(0), 

EDI fixed statie init(0), 

CF float ext, 

AE{8) float ext, 

DCLOCK fixed ext, 
eode fixed(7) , 
result ptr, 
m ( 5) label , 

( ia_ptr , ib_ptr ,ie ptr , id_ptr , ie_ptr ) ptr; 
^inelude 'del inp . txt"^*» 

%inelude 'input. ttt'; 
if eode = 1 then do; 

if EE < ECLOCK then do ; 

if (nr_inputs> 0) & (nr_inputs< 6) 
then do; 
m { 1 ) =la ; 
m(2)=lb; 
m (3 )=lc; 
m(4)=ld ; 
m (5)=le; 

goto m( nr_inputs ) ; 
end; then / 
else goto If; 

le: ie_ptr = input (5) (eod(5)); 

Id: id_ptr = input (4) (cod(4)); 

le: ic_ptr = input (3) (cod(3)); 

lb: ib_ptr = input (2) (eod(2)); 

la: ia_ptr = input (1) (cod(l)); 

If: 0 = ?; 

EEl = ECLOCK; 

%inelude 'fune.txt'; 

• 

EE = ECLOCK; 
end; /=5= if EE -/ 
end; /- if code then 
else do; 

if EEl < ECLOCK then do; 

0 = f; 

EEl = ECLOCK; 
end; /* EEl */ 
end; else */ 

:iinclude 'out. ttt'; 
result = addr(F); 
return (result); 
end ; 
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/* DA.PLI ’^/ 

^include 'dclfun . ttt 



del 

ftc entry (float binary) returns (char (17 ) var); 

del 

TIME fixed ext, 

T float ext, 

T1 float ext, 

input(5) variable entry ( fixed '7) ) r eturns ( ptr ) , 
CF float ext, 

AD(8) float ext, 

EA(8) float ext, 
code f ixed(7) , 
result ptr, 
m ( 5 ) label , 

(ia_ptr,ib_ptr,ic_ptr,id_ptr,ie_ptr) ptrj 

^include "del inp . txt ' J 
%include "input. ttt"; 



if TIME > 0 then ao ; 

if (nr_inputs> 0) & (nr_inputs< 6) 
then'do »* 
m ( 1 ) =la ; 
m(2)=lb; 
m (3 )=lc ; 
m(4)=ld; 
m( 5 )=le ; 

i?oto m(nr inputs ) ; 
end; /* then */ 
else goto If; 



le ; 


ie_ptr 


= input 


(5) 


(1); 


Id: 


id_ptr 


= input 


(4) 


(1); 


Ic: 


ic_ptr 


= input 


(3) 


(1); 


lb: 


ib_ptr 


= input 


(2) 


(1); 


la : 


ia~ptr 


= input 


(1) 


(1); 


If : 










%include 


"f unc . tx t " ; 


/- 



end; /* then 
result = addr(F); 



IB - ECT */ 



)J,include "out.ttt"; 
J 

return (result); 
end ; 
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REFORM: proc opt i ons (ma in ) »' 
del 

( input , output ) file, 
buff char(l) varying; 

open file(input) stream env(b(1024)) title( '$1 .$1 ') ; 

open file(output) stream output env(b(1024)) title( '$2 .$2') 

do while ( "l 'b ) * 

read file (input) into (buff); 
buff = translate (buff ; 
write file (output) from (buff); 
end ; 

end reform; 



PRINT: proc options (main) ; 
del 

da ta dec (15,6) , 

( j ,k,m) fixed bin (7) , 

( time ,dat, temp_str ) char(128) var, 
( inf 0 , input , output ) file, 

^include 'datainfo.txt'; 



open file(input) record sequential t i t le ( 'EAT 1 . FLT ' ) 
env ( b (512 ) ) ; 

open file(output) stream output ti tie ( '$1 1 ' ) 

env(b(512) ) ; 

do k = 0 to timei; 
time = ' '; 

temp_s tr = 'T IME= ' 1 'char(k); 
do m = 1 to length( temp_str ) ; 

if substr ( temp_str ,m, 1) ~= ' ' then 
time = time 1 1 substr( temp_str ,m ,1 ) ; 

end; 

time = '"'j ^m 'lit ime ; 

write file(output) from(time); 
do j= 1 to nr_data; 

read file^input) into(data); 

dat = ' 1 1 name ( ,i ) 1 1 '= ' 1 1 cha r (da ta ) ; 

write file(output) from(dat); 
end ; 

end; /•'* do k 
end ; 
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aoilindHV * * * ' ' aaniiidHV 



ArPENIlX C 



GRAPHS OF EXAy;PL£ PROBLEM RESULTS 




Figure ISA 

Response to Unit Step Input With No Correction 




+ 



Figure IS-B 

Response to Unit Step Input *'ith Correction 
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amplitude AMPLITUDE 




TIME 



Figure 19C 

Terminated Bamp Input With No Corrections 




Figure 19D 

Terminated Ramp Input With Corrections 
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AMPLITUDE AMPLITU DE 





rigure ISE 

Response to Variable Step Inputs lith and Without Corrections 



AMPLITUDE amPLI 





figure l&F 

Response to Variable Step Inputs *itb and Without Corrections 
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